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SECTION  1 


INTRODUCTION 


This  report  is  a  comprehensive  user's  guide  to  the  CM3DT  blunt 
body  code  for  the  calculation  of  inviscld,  hypersonic  flow  over  ablated 
reentry  vehicle  nosetips.  This  code  provides  numerical  solutions  to  the 
exact  Inviscld  time -dependent  fluid  mechanic  equations  and  is  formulated 
in  terms  of  a  generalized  coordinate  system  that  is  automatically  aligned 
with  asymmetric  nosetlp  geometries.  The  analysis  from  which  this  code  has 
been  developed  may  be  found  in  Reference  1. 

A  brief  description  of  the  techniques  used  in  this  code  is 
given  in  Section  2.0,  and  the  logical  organization  of  the  code  is  described 
In  Section  3.0.  The  inputs  required  by  this  code  are  defined  in  Section 
4.0,  and  Section  5.0  describes  the  output  generated. 

The  coupling  of  this  nosetlp  flow  field  code  to  an  inviscid 
supersonic  afterbody  code,  enabling  inviscid  aerodynamic  predictions  to 
be  made  for  entire  reentry  vehicles,  is  described  in  Section  6.0.  Samples 
of  the  code  Inputs  are  provided  In  Section  7.0.  Section  8.0  presents 
criteria  to  define  those  Indented  nosetlp  shapes  for  which  the  CM3DT  code 
can  be  expected  tu  produce  accurate  inviscid  flow  field  results. 

The  CM3DT  code  was  developed  on  the  CDC  Cyber  176  computer,  and 
is  compatible  with  the  CDC  7600  computer.  Qualified  users  may  obtain  this 
program  through  BMO/SYDT. 
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SECTION  2 

DESCRIPTION  OF  THE  TECHNIQUE 


The  CM3DT  code  is  a  numerical  procedure  for  the  calculation 
of  hypersonic  invlscid  flows  over  reentry  vehicle  nosetlps.  This 
procedure  has  been  developed  with  the  goal  of  computing  flows  over  ablated, 
asyntr«tr1c  nosetlp  geometries  which  could  not  be  treated  by  previous 
Invlscid  blunt  body  codes.  In  particular,  this  technique  has  extended 
the  range  of  nosetlp  geometries  for  which  successful  Invlscid  flow  field 
calculations  can  be  performed  by  using  a  new  coordinate  system  which  can 
be  closely  aligned  with  the  body  geometry  and  by  using  the  x-dlfferencing 
scheme  for  calculations  on  Indented  shapes  where  embedded  shocks  form. 

The  CM3DT  code  computes  the  steady  flow  over  a  nosetlp  as  the 
asymptotic  limit  of  an  unsteady  flow.  Integrating  the  Invlscid  time- 
dependent  equations  with  a  second-order  accurate,  forward-marching  (in 
time),  explicit  finite  difference  scheme.  Body  and  bow  shock  points  are 
computed  using  forms  of  the  appropriate  characteristic  compatibility  condi¬ 
tions.  Because  of  the  Invlscid  flow  assumption,  this  technique  is  limited 
to  high  Reynolds  number  flows  (thin  boundary  layers)  without  separation. 

'•lie  coordinate  system  used  in  this  technique  is  the  result  of 
conformal  transformations  carried  out  in  each  meridional  plane,  closely 
aligning  one  of  the  coordinates  with  the  body,  and  having  another  coordinate 
nearly  normal  to  the  body.  The  conformal  transformations  are  defined  in 
terms  of  "hinge  points",  which  are  discrete  points  selected  within  each 
meridional  plane  that  approximate  the  body  meridional  curve.  The  hinge 
points  may  be  specified  by  the  user  of  the  technique,  or  can  be  generated 
internally  by  the  code. 
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Numerical  calculations  are  performed  In  a  computational  coordi¬ 
nate  system.  In  which  mesh  points  are  equally  spaced  In  the  transformed 
coordinate  system:  circumferentially,  along  the  body  profile,  and  between 
the  body  and  the  shock. 

Two  computational  procedures  are  available  within  the  CM3DT 
code;  the  basic  non-conservation  approach  and  the  X-dlfferencing  approach. 
The  non-conservation  procedure  Is  the  standard  technique  for  most  convex 
nosetlp  flow  field  calculations;  however,  for  calculations  where  embedded 
shocks  form,  the  X-dlfferenclng  technique  must  be  used.  The  X-dlfferenclng 
technique  Is  capable  of  computing  embedded  shocks  In  an  approximate  manner 
by  use  of  a  finite  difference  technique  that  accurately  models  the  physical 
domain  of  dependence  of  each  point  being  computed. 

The  CM3DT  blunt  body  code  has  been  coupled  to  a  steady  super¬ 
sonic  afterbody  code  to  enable  the  prediction  of  invlscid  aerodynamic 
characteristics  for  an  entire  reentry  vehicle  with  an  ablated  nosetlp. 

For  a  maximum  flexibility,  the  nosetlp  centerline  may  be  rotated  In  both 
pitch  and  yaw  relative  to  the  afterbody  axis,  which  may  Itself  be  bent  in 
the  pitch  plane. 

For  a  typical  blunt  body  calculation,  with  6  points  between 
the  body  and  bow  shock,  18  points  along  the  body,  and  9  circumferential 
planes,  and  run  for  400  time  steps,  the  non-conservation  code  requires 
approximately  two  minutes  central  processor  time  on  the  CDC  Cyber  176 
computer,  while  the  x-differenclng  procedure  requires  approximately  three 
minutes.  (In  comparison,  the  blunt  body  code  formulated  in  a  spherical 
coordinate  system  described  in  References  2  and  3  requires  slightly  more 
than  one  and  a  half  minutes  for  the  same  esse.)  The  CM3DT  code  requires 
81000  words  of  core  storage. 

The  analysis  on  which  the  CM3DT  coda  Is  based  is  described 
In  more  detail  In  Reference  1. 
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SECTION  3 


LOGICAL  ORGANIZATION  OF  THE  CODE 


This  section  contains  a  description  of  the  logical  flow  of  the 
CM3DT  program.  An  alphabetical  listing  of  the  subroutines  found  in  CM3DT, 
together  with  a  brief  description  of  each  routine,  may  be  found  in  the 
Appendix,  Section  10.3. 

The  main  program  (CM3DT)  reads  and  prints  out  the  Input  data. 
Freestream  and  stagnation  point  conditions  are  calculated  (FREE,  ATMP, 
STAGPT).  For  a  real  gas  calculation,  tables  of  shock  properties  are 
computed  (SHKTAB) .  Hinge  point  locations  are  calculated,  if  not  input, 
and  printed  out  (HING3D) ,  The  required  geometric  parameters  are  deter¬ 
mined  using  GE0M3  (analytic  geometry),  GE0M4  (arbitrary  geometry  with  no 
centerline  offset),  or  GEQM5  (arbitrary  geometry  with  centerline  offset). 
An  approximate  bow  shock  shape  is  defined,  the  computational  grid  is 
determined,  and  coordinate  transformation  parameters  are  calculated  for 
each  grid  point  (GRID3D,  MAP3D).  Initial  estimates  of  pressure,  entropy 
and  the  transformed  velocity  components  are  made  at  each  point  In  the 
computational  grid.  If  the  option  to  save  data  on  a  binary  file  has  been 
chosen,  the  binary  file  is  initialized  and  thQ  data  written  (BFL3CM). 

(In  the  case  of  a  restart,  all  Initial  data  are  read  from  the  binary 
file.) 


The  difference  equations  for  the  flow  variables  of  pressure, 
velocity  components,  and  entropy  are  solved  in  PVE3DNC  If  the  non¬ 
conservation  form  has  been  chosen,  or  in  PVE3DL  If  the  X-scheme  is 
desired.  At  the  ginning  of  each  of  these  subroutines,  the  time  step 
counter  Is  incremented  by  one.  When  the  predictor-corrector  loop  for  a 
time  step  has  been  completed,  data  at  each  grid  point  are  printed  if 
output  is  desired  at  that  step,  and  saved  on  the  binary  file  if  at  a  step 
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specified.  If  the  current  time  step  Is  not  the  last  step  desired,  a 
transfer  back  to  the  beginning  of  PVE3DNC  or  PVE3DL  Is  irudej  otherwise, 
control  returns  to  CM30T,  and  the  program  terminates. 


The  CM3DT  program  requires  81000  (decimal)  words  of  memory  for 
execution,  without  overlays.  The  program  may  be  compiled  and  executed 
from  an  UPDATE  program  library  or  loaded  and  executed  from  a  user-created 
LIBRARY.  In  the  first  case*  It  Is  necessary  to  omit  from  the  list  of 
subroutines  to  be  compiled  either  PVE3DNC  (If  the  X-scheme  option  Is 
chosen),  or  PVESu.  (If  the  non-conservation  method  Is  desired).  Since 
these  subroutines  are  mutually  exclusive,  loading  both  would  unnecessarily 
Increase  the  core  storage  required.  In  like  manner,  loading  from  a  LIBRARY 
requires  use  of  the  LDSET(0MIT-PVE3DNC)  or  LDSET(QMIT*PVE3DL )  comnand  In 
the  CDC  control  card  stream.  Regardless  of  the  method  of  loading  the 
program,  It  is  necessary  to  preset  core  memory  to  zero  before  Initiating 
execution.  Following  are  examples  of  each  of  the  above  methods  of  loading 
and  executing  CM3DT: 


1.)  Compiling  and  executing  CM3DT  from  an  UPDATE  program  library: 


JOBCARD? 

ACCOUNT  * 

REQUEST ?  TAPE21 *  ♦PF. 

ATTACH? OLBPL  ?  CM3DTPL? I D« _ 

UPDATE. 

ftn  a -compile? b»cm3dtb?l»o:> 
ldsetlpreset*zero> 

CM3DTB. 

CATALOG?  TAPE21  ?  CM3DTDATA?  RP-9991?  ID- _ _ _ 

7' 8 

♦COMPILE  CM3DT ?  PVE3DNC?  HING3D?  GRIDS!'? 0UT3B?  BFL3CM?  MAP3D 
♦COMPILE  SfcQMS? GE0M4? GE0M5? STAGPT ? SHKTAB?  SHOCK?  FREE? RGAS? ATMP 
♦COMP I LE  CUFT 1 • CUFT£ ?  TBL1 ? ATERR? I  DEAL ?  L I HB  ?  RLFPR  ?  TBL ?  TLU 1 

7.--S--9 

CM3DT  INPUT  DATA 

6s?s8s9 
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2.)  Loading  and  executing  CM3DT  from  a  user  library: 


ACCOUNT, 

3TT$CH* L  IB*  CMSBTLIBs  ID« _ 

LIBIWWY<LIB> 

LBiET  <»RESE  WE*0» On  I T«PVE3DNC > 
LIBLQ3B  «H  B » Cri'SDT  > 

EXECUTE. 

?  -  5  •? 

CMSBT  IN»UT  DPlT* 

S' 9 


The  file  specified  for  storage  of  binary  Information  In  the  CM30T 
program  Is  TAPE21.  Since  TAPE21  Is  both  an  Input  and  an  output  device  to 
the  CM3DT  program,  It  is  advisable,  in  cases  of  restart,  for  the  user  to 
attach  the  data  file  from  which  the  program  Is  to  be  restarted,  with  some 
designator  other  than  TAPE21,  and  then  copy  this  file  to  a  TAPE21  local 
to  the  job,  to  eliminate  loss  of  data  in  the  event  of  error. 

In  coupling  the  CM3DT  program  to  the  afterbody  code,  it  is  first 
necessary  to  perform  the  CM3DT  calculation  and  save  the  data.  Then  the 
afterbody  program  must  be  modified  In  the  following  manner: 

1. )  Subroutines  CM3ST,  REALP,  and  BFL3CH  must  be  in¬ 

corporated  into  tho  afterbody  code, 

2. )  The  variables  XO,  ZO,  ZB  and  TE$T(4)  must  be  added 

te  the  list  of  Input  variables,  with  the  default 
values  of  the  first  three  set  at  0.0  and  the 
fourth  at  I.0E-3. 
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SECTION  4 


CODE  INPUTS 


This  section  details  the  Input  variables  required  by  the 
CM3DT  code.  The  Input  variables  are  discussed  In  thel/  natural  groupings; 
l.e.,  geometry  specification,  freestream  conditions,  grid  selection,  out¬ 
put  controls,  and  special  options.  A  glossary  of  all  Input  variables  Is 
provided  In  the  Appendix,  Section  10. I. 

4.1  GEOMETRY 

The  body  geometry  can  be  specified  In  either  of  two  ways: 
analytically  (for  axlsymmetrlc  geometries)  or  by  tabulation  of  body  points 
(for  arbitrary  geometries).  Several  options  are  available  within  the 
arbitrary  geometry  definition  procedure  to  simplify  the  Input  required  for 
special  cases;  these  options  are  detailed  In  Section  4.1.2. 

Several  geometric  Input  parameters  are  common  to  both  the 
analytic  and  arbitrary  geometric  definition  procedures.  In  both  cases 
the  geometry  Is  defined  relative  to  a  user  defined  nosetlp  geometric 
axis,  which  must  Intersect  the  afterbody  axis.  For  generality,  the  after¬ 
body  axis  may  be  bent  In  the  pitch  plane,  and  the  nosetlp  geometric  axis 
may  be  rotated  In  both  pitch  and  yaw  from  the  bent  afterbody  axis. 

The  orientation  of  the  nosetlp  geometric  axis  relative  to  the 
afterbody  axis  Is  Illustrated  In  Figure  4.1.  The  bend  angle  (If  any)  of 
the  afterbody  axis  In  the  pitch  plane  Is  denoted  by  6^  (DELA),  Input  In 
degrees.  The  angles  6^  (DELI)  and  <Sg  (DEL2) ,  also  specified  In  degrees, 
define  the  orientation  of  the  nosetlp  geometric  axis  relative  to  the  bent 
afterbody  axis  In  the  pitch  and  yaw  planes,  respectively.  This  procedure 
for  the  definition  of  the  nosetlp  geometric  axis  (x-axIs)  allows  the  user 
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the  flexibility  of  using  the  most  appropriate  axis  for  defining  an  ablated 
nosetlp  geometry,  Independent  of  the  orientation  of  the  nosetlp  shape 
relative  to  the  afterbody.  The  geometric  nosetlp  axis  Intersects  the  bent 
afterbody  axis  at  x  8  xQ(XH0). 

Two  other  geometric  parameters  must  be  specified  in  all  cases. 

The  first,  XBO,  is  the  point  on  the  x-axis  that  intersects  the  body  surface, 
as  shown  in  Figure  4.2.  The  second  parameter,  the  airay  XBB(L),  specifies 
the  x  locations  at  the  body  that  define  the  downstream  boundary  of  the 
region  to  be  computed,  as  Illustrated  In  Figure  4.2.  Details  on  the  use 
of  the  XBB  array  vary  with  the  geometry  option  used,  and  are  discussed  In 
Sections  4.1.1  and  4.1.2. 

In  all  cases,  the  values  in  the  XBB  array  must  be  selected 
sufficiently  far  downstream  that  the  flow  across  the  downstream  boundary 
is  supersonic.  Additionally,  the  region  of  Influence  of  any  downstream 
boundary  point  must  not  contain  any  other.  The  downstream  boundary  must 
be  selected  so  that  the  region  computed  will  include  complete  shock  layer 
information  In  the  initial  plane  for  the  afterbody  calculation. 

4.1.1  Analytic  Geometry 

Sufficiently  regular,  axisymmetric  geometries  may  be  defined 
analytically  in  CM3DT.  The  analytic  description  consists  of  specification 
of  coefficients  of  a  general  equation  for  each  of  several  curve  segments. 
Wi*h  this  geometry  option,  the  nosetlp  geometric  coordinate  system  (x,y,$) 
serves  as  the  computational  coordinate  system  (x,y,4>).  The  number  of 
segments  to  be  used  is  specified  by  NSEG;  a  maximum  of  six  segments  is 
allowed.  Endpoints  of  the  segments  are  defined  in  the  ADIV(J)  array; 

NSEG-1  values  of  ADIV  must  be  input.  Within  each  segment,  the  curve  is 
defined  from 

y(x)  s  A:  +  A2x  +  A3x2  +  A4x3  +  AgX4  +  AgX5 

+  A?x  Aa  +  A9  /A102  -  (x  -  An)2  (4.1) 
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where  x  *  x  -  A^.  The  coefficients  for  each  segment  j  are  input  via 
the  AC0F(I,J)  array. 

When  using  the  analytic  geometry  option,  the  user  must  ensure 
that  the  analytic  description  Is  consistent  with  the  value  of  X30  specified, 
and  that  the  segments  specified  are  continuous. 

To  ceflne  the  downstream  boundary  of  the  computation  region  when 
using  this  geometry  option,  only  XBB(l)  need  be.  input,  as  the  downstream 
boundary  is  assumed  to  be  axi symmetric  for  axi symmetric  nosetlp  geometries. 

,  i  ■ 

4.1.2  Arbitrary  Geometry 

Two  arbitrary  geometry  options  are  available  within  the  CM3DT 
code.  In  both  cases,  the  nosetlp  geometry  is  defined  In  the  (x,y,$) 
cylindrical  coordinate  system,  which  intersects  the  bent  afterbody  axis 
at  x  *  xQ  (XHO)  and  is  oriented  relative  to  the  bent  afterbody  axis  by  the 
angles  (DELI)  and  dg  (DEL2).  The  first  arbitrary  geometry  option  allows 
definition  of  profiles  in  meridional  planes  to  establish  the  nosetlp  geometry 
(with  several  special  options  available  for  defining  the  nature  of  the 
cross-sections).  The  second  option  allows  the  specification  of  cross- 
sections  (of  several  special  types)  at  specified  axial  stations,  with  the 
cross-section  centers  potentially  offset  from  the  x  axis.  These  options  are 
described  in  greater  detail  below. 

4,1.2. 1  Arbitrary  Geometry  Without  Centerline  Offsets  (NPHI>0) 

With  this  option,  arbitrary  geometries  are  defined  by  tabular 
specification  of  points  XB(I,L),  YB(I,L)  in  the  $  planes  defined  by  the 
array  PHIBD(L)  (degrees).  NPHI  (which  is  a  positive  number  for  this  option) 
specifies  the  number  of  PHIBD  planes  in  which  XB,  YB  points  are  to  be 
defined,  and  the  array  NPTS(L)  specifies  the  number  of  XB,  YB  points  input 
in  each  plane.  A  maximum  of  nine  planes  of  geometry  may  be  defined,  and  a 
maximum  of  50  XB,  YB  points  may  be  specified  in  each  plane. 
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When  this  arbitrary  geometry  option  Is  used,  the  computational 
(x,y,<j>)  coordinate  system  Is  Identical  to  the  nosetlp  geometry  (x,y,$) 
coordinate  system.  The  program  automatically  sets  XB(1,L)  to  XBO  and 
YB(1,L)  to  0. 

Cubic  spline  fits  of  the  XB,YB  data  are  used  both  longitudinally 
and  circumferentially  to  determine  the  body  geometry  at  any  point.  A 
number  of  options  are  available  with  the  arbitrary  geometry  procedure  to 
provide  for  special  cases  In  which  the  cross-sectional  character  of  the 
nosetlp  is  known.  These  cases  are  enumerated  below,  and  are  Illustrated 
in  Figure  4.3.  Use  of  these  special  cases  is  determined  by  the  value  of 
NPHI  specified. 


This  option  Is  exercised  by  setting  NPHI  *  1.  Circular  cross- 
sections  are  assumed,  and  only  one  XB.YB  table  need  be  input,  along  with 
NPTS(l).  The  PHIBD  array  is  not  required  for  this  ca*'®. 


Offset  Circular  Cross-Sections  (NPHI  *  2) 


With  this  option,  the  body  cross-sections  are  assumed  to  be 
circular,  with  the  centers  of  the  cross-sections  offset  from  the  nosetlp 
computational  axis  in  the  pitch  plane.  Use  of  this  option  requires  setting 
NPHI  =  2  and  supplying  two  sets  of  XB.YB  data,  with  the  appropriate  values 
of  NPTS.  The  PHIBD  array  Is  automatically  set  to  0°,  180°,  and  no  PHIBD 
Input  is  required  of  the  user.  The  arrays  XB ( I , 1 )  and  YB(I,1)  are  assumed 
to  define  points  In  the  $  *  0°  plane,  and  the  XB(I,2)  and  YB(I,2)  arrays 
are  assumed  to  define  points  In  the  $  »  180°  plane,  relative  to  the  nosetlp 
geometric  coordinate  system. 
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Bl-Elllptlc  Cross-Sections  (NPHI  *  3) 


By  setting  NPHI  *  3,  bi-ell iptlc  cross-sections  are  assumed, 
and  the  body  geometry  Is  assumed  to  be  symmetric  about  the  pitch  plane. 
Three  planes  of  XB.YB  data  with  their  associated  values  of  NPTS  are  re¬ 
quired  and  are  assumed  to  apply  to  the  $  *  0°,  90°,  and  180°  planes. 

With  this  option,  each  quadrant  of  the  cross-section  Is  fit  as  the 
quadrant  of  an  ellipse.  Again,  the  PHIBD  array  Is  not  required  for  this 
option. 


Quad-Elliptic  Cross-Sections  (NPHI  *  4  and  ISA  ■  2) 

The  quad-elliptic  cross-section  option  is  exercised  by  setting 
NPHI  =  4  with  ISA  »  2  (no  pitch  plane  of  symmetry).  With  this  option, 

XB, YB  and  NPTS  data  must  be  input  in  the  <J>  *  0°,  90°,  180°,  and  270° 
planes  and  no  PHIBD  input  Is  required.  Each  quadrant  of  the  cross-section 
is  assumed  to  be  the  quadrant  of  an  ellipse.  (Note  that  when  NPHI  *  4 
and  ISA  =  1,  arbitrary  cross-sections  with  a  pitch  plane  of  symmetry  are 
assumed. ) 


When  using  this  arbitrary  geometry  option  (NPHI>0),  the  down¬ 
stream  boundary  of  the  computational  region  Is  defined  by  specifying  the 
array  XBB(L) ,  where  each  element  defines  the  location  of  the  downstream 
boundary  at  the  body  In  the  corresponding  PHIBD(L)  plane  of  input  geometry. 
The  appropriate  values  of  XBB  In  the  planes  to  be  computed  are  determined 
automatically  within  CM3DT  by  linear  interpolation  in  $. 
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4. 1.2. 2  Arbitrary  Geometry  with  Centerline  Offsets  (NPHI<0) 


In  using  this  arbitrary  geometry  option,  body  cross-sections  of 
certain  limited  classes  are  defined  In  specified  planes  normal  to  the 
nosetlp  geometric  axis,  x.  The  center  of  each  cross-section  may,  however, 
be  offset  In  both  the  pitch  and  yaw  planes  from  the  x  axis.* 

This  geometry  option  Is  Invoked  by  specifying  a  value  of  NPHI<0; 
the  value  of  NPHI  determines  the  type  of  cross-section  to  be  defined,  as 
shown  In  Figure  4.4.  Admissible  cross-sections  Include  circular,  elliptic, 
bi -elliptic,  and  quad-elliptic. 

Stations  at  which  the  cross-sections  are  to  be  defined  are 
specified  through  the  array  XE(I),  where  XE  represents  x,  the  axial  co¬ 
ordinate  In  the  nosetip  geometric  coordinate  system  (which  is  related  to 
the  bent  afterbody  axis  through  the  values  input  for  DELI,  DEL2,  and  XHO). 
Cross-sections  may  be  defined  at  up  to  50  stations,  and  XE(l)  must  equal 
XBO.  The  number  of  cross-sections  input  is  set  by  NPTS(l). 

At  each  of  the  specified  stations,  the  offset  of  the  cross- 
section  center  from  the  geometric  axis  is  defined  by  the  arrays  Dl(I) 
(offset  in  pitch  plane)  and  D2(I)  (offset  in  yaw  plane),  as  shown  in 
Figure  4.5.  An  option  also  exists  in  which  the  centerline  offset  in  the 
pitch  plane  may  be  defined  analytically  through  the  ACOF  and  ADIV  arrays, 
using  the  same  approach  as  for  the  analytic  description  of  axi symmetric 
body  geometries,  except  that  the  number  of  segments  used  is  input  as  a 
negative  number  (NSEG<0).  With  the  pitch  offset  defined  analytically,  the 
yaw  offset  is  assumed  to  be  0. 

Parameters  describing  each  cross-section  are  input  through  the 
arrays  AE'I),  BE(I),  CE(I),  and  DE(I),  which  correspond  to  the  quantities 
a,  b,  c,  and  d  shown  in  Figure  4.4.  The  types  of  crrss-sectlons  allowed 
are  detailed  below. 

*Nott  that  if  <5,  0,  f  0,  or  6.  i  0,  the  nosetlp  cross-sections  will  not 

be  p?rallel  to  afterbody  cross-sections.  It  is  the  responsibility  of  the 
user  to  ensure  compatibility  between  the  nosetip  and  afterbody  geometry 
definitions. 
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Circular  Cross-Sections  (NPHI  •  -1) 


This  option  Is  exercised  by  setting  NPHI  ■  -1.  Circular  cross- 
sections  are  assumed,  and  Input  Is  required  for  the  XE  and  AE  arrays  only, 
as  well  as  the  definition  of  the  offset  of  each  cross-section  center. 

Elliptic  Cross-Sections  (NPHI  »  -2) 

By  setting  NPHI  ■  -2,  elliptic  cross-sections  are  assumed,  with 
Input  required  for  the  XE,  AE,  and  BE  arrays,  plus  definition  of  the  offset 
of  the  cross-section  centers. 

Bl-Elllptlc  Cross-Sections  (NPHI  ■  -3) 

With  this  option,  the  body  cross-sections  are  assumed  to  be  bi- 
elllptic.  Use  of  this  option  requires  setting  NPHI  »  -3  and  supplying 
values  for  the  arrays  XE,  AE,  BE,  and  CE,  as  well  as  the  Information 
required  to  de.lne  the  offset  of  the  cross-section  centers. 

Quad-Elliptic  Cross-Sections  (NPHI  ■  -4) 

The  quad-elliptic  cross-section  option  os  exercised  by  setting 
NPHI  «  -4  and  supplying  data  for  the  XE,  AE,  BE,  CE,  and  DE  arrays,  as 
well  as  the  definition  of  the  offset  of  each  cross-section  center. 

• 


When  this  arbitrary  geometry  option  is  used,  it  is  necessary 
to  define  a  computational  (x,y,$)  coordinate  system  that  differs  from  the 
geometric  (x,y,$)  system  to  avoid  multi-valued  body  geometric  functions 
y  ■*  yb(x,$).  The  computational  x-axIs  Is  defined  automatically  within 
CM3DT  as  a  straight  line  that  passes  through  the  most  forward  point  of 
the  body,  specified  by  XE(1)  *  XBO,  01(1),  and  D2(l),  and  that  intersects 
the  x  axis  at  xQ  (XHO).  (Note  that  the  Input  XHO  in  CM3DT  is  identical 


to  the  Input  XO  In  the  Afterbody  code,  described  in  Section  7.2.)  Details 
of  the  coordinate  transformations  are  provided  In  the  Appendix,  Section 
10.4. 


Calculations  required  for  the  computational  centerline  are 
carried  out  Internally  within  CM3DT,  Including  redefinition  of  XBO  and 
XHO  for  consistency  with  the  new  coordinate  system,  with  XBO  set  equal  to 
zero  at  the  geometric  stagnation  point.  In  addition,  new  values  of 
and  62  are  computed  to  define  the  orientation  of  the  new  computational  axis 
relative  to  the  bent  afterbody  axis.  New  values  of  a  and  B  (relative  to 
the  computational  axis)  are  also  computed  from  the  new  values  of  6j  and 
o2,  as  discussed  In  Section  6. 

When  using  this  geometry  option,  the  downstream  boundary  of  the 
computational  region  Is  defined  by  specifying  XBB(l);  the  downstream 
boundary  Is  then  taken  to  be  the  axial  location  XE(I)  that  Is  the  first 
one  located  downstream  of  XBB(l).  This  boundary  In  the  x-system  Is  then 
transformed  Internally  within  the  CM3DT  code  to  produce  a  new  array  of 
XBB  data  for  the  planes  to  be  computed  In  the  computational  (x,y,4>) 
coordinate  system. 

4,2  FREESTREAM  CONDITIONS  AND  THERMODYNAMICS 

4.2.1  Freestream  Conditions 

For  each  calculation  It  Is  necessary  to  specify  the  freestream 
Mach  number  (AMINF)  or  the  freestream  velocity  (QINF),  as  well  as  the 
static  pressure  and  temperature  (PINF  and  TINF).  Values  of  PINF  and  TINF 
will  be  generated  automatically  from  the  1962  standard  atmosphere  tables 
by  specifying  the  desired  altitude  (ALT)  and  the  flag  IATMP  f  0.  For 
non-dimensional  Ideal  gas  calculations,  PINF  and  TINF  are  set  to  1.0 
(default  values). 
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The  orientation  of  the  vehicle  relative  to  the  freestream 
velocity  vector  Is  defined  by  specification  of  the  angle  of  attack  (ALPHA) 
and  the  sldesl Ip  angle  (BETA).  The  angles  ALPHA  and  BETA,  specified  In 
degrees,  are  defined  relative  to  the  unbent  afterbody  axis:  the  CM3DT 
code  will  automatically  compute  the  values  of  a  and  0  relative  to  the  nose 
tip  computational  axis  as  required  for  the  computational  procedure.  Thus 
the  orientation  of  the  nosetlp  centerline  relative  to  the  wind  vector  Is 
defined  In  terms  of  ALPHA,  BETA,  DELA,  DELI,  and  DEL2.  (Note  that  If 
NPMI<0,  the  values  of  DELI  and  DEL2  are  recomputed  Internally  by  CM3DT.) 

When  ISA  •  1  Is  Input,  signifying  a  pitch  plane  of  symmetry, 
the  sideslip  angle  BETA  Is  automatically  set  to  zero  by  the  code. 

Sign  conventions  for  ALPHA  and  BETA  are  shown  In  Figure  4.6; 
positive  ALPHA  Implies  nose  up  and  positive  BETA  Implies  nose  left, 
relative  to  the  freestream  velocity  vector. 

4.2.2  Thermodynamics 

Two  options  are  available  within  the  CM3DT  code  for  the 
thermodynamic  model  used:  Ideal  gas  and  equilibrium  air  (real  gas). 

These  options  are  described  In  more  detail  below. 

4. 2. 2.1  Ideal  Gas  (IRG  »  0) 

The  standard  operating  procedure  for  Ideal  gas  calculations  Is 
to  use  non-dimensional  variables,  In  which  PINF  *  TINF  *  RIDEAL  ■  1. 

With  this  normalization  the  variables 

CM3DT 

P 

P 

T 

(u.v.w) 
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take  the  forms: 

ACTUAL 

P/P. 

P/P» 

T/T. 

(u.v.w)/^ 


W.« 


Note  that  with  this  formulation  the  gas  law  becomes  p  *  pT.  The  isentropic 
exponent  is  specified  as  GIDEAL.  If  desired,  a  compressibility  factor 
(ZIDEAL)  may  also  be  specified. 

Dimensional  ideal  gas  caluclations  are  also  possible,  requiring 
only  that  a  dimensionally  consistent  value  of  RIDEAL  Is  specified,  as  well 
as  appropriate  values  of  PINF  and  TINF. 

"'1. 2.2.2  Real  Gas  (IRG  =  1) 

/ 

/ 

All  quantities  in  real  gas  calculations  are  dimensional  with 

the  units: 

CM3DT  ACTUAL 

p  lbf/ft2 

p  slugs/ft3 

T  *R 

(u,v,w)  ft/sec 

The  real  gas  thermodynamic  properties  available  in  CM3DT  are  for 
equilibrium  air.  Gas  properties  required  during  the  calculation  are  ob¬ 
tained  by  interpolation  on  the  thermodynamic  property  tables  supplied  with 
the  code. 

4.3  GRID  DEFINITION 

In  the  CM3DT  code  a  conformal  transformation  is  used  to  map  the 
cylindrical  (x ,y ,0)  computational  space  onto  (S»n,9)  space,  where  n  *  con¬ 
stant  surfaces  are  closely  aligned  with  the  bpdy  surface  and  5  *  constant 
surfaces  are  nearly  normal  to  the  surface.  The  transformation  used  is 
based  on  a  number  of  "hinge"  points  selected  within  each  $  plane  being 
computed  that  lie  within  the  body  profile  and  closely  model  the  body  shape. 
These  hinge  points  may  either  be  input  by  the  user  or  are  generated 
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automatically  by  the  code,  except  when  the  offset  centerline  geometric 
option  Is  used  (In  which  case  the  hinge  points  will  be  generated  auto¬ 
matically  In  the  computational  coordinate  system).  The  conformal  trans¬ 
formation  will,  In  general,  be  non-ax  1  symmetric  (IMAP  •  1;  default  value), 
unless  the  user  specifies  an  axisymmetric  transformation  (IMAP  *  0)  or 
the  nosetlp  geometry  Is  axisynmetric. 

4.3.1  Automatic  Generation  of  Hinge  Points 

The  standard  option  In  CM3DT  is  to  have  the  hinge  points  generated 
automatically,  selected  by  setting  IHINGE  ■  0  (default  value).  The  hinge 
points  to  be  generated.  Illustrated  In  Figure  4.7,  consist  of  one  point 
on  the  centerline  in  front  of  the  body,  one  on  the  centerline  inside  the 
body,  and  the  rest  equally  spaced  in  body  wetted  length,  located  inside 
the  body  along  inward  body  normals.  The  last  hinge  point  is  located 
downstream  of  the  end  of  the  body  geometry  to  ensure  a  smooth  transformation. 
If  the  offset  centerline  geometric  option  is  used,  the  hinge  points  must  be 
generated  automatically. 

The  number  of  hinge  points  generated  in  each  plane  is  determined 
by  the  input  JA;  JA+2  hinge  points  will  be  generated  in  each  plane.  The 
maximum  value  of  JA  is  7.  If  an  axisymmetric  body  geometry  is  specified, 
the  hinge  point  locations  will  be  the  same  in  each  <j>  plane  (axisymmetric 
mapping) . 


The  first  hinge  point  is  located  at  a  distance  of  TEST  (8)  *  Afl 
forward  of  the  body,  where  AQ  is  an  estimate  of  the  shock  stand-off  distance 
computed  automatically  within  the  code  as  a  function  of  Mw  and  RN,  the 
effective  nose  radius.  Subsequent  hinge  points  are  located  a  distance 
of  DELH*RN  along  inward  body  normals  equally  spaced  in  wetted  length  along 
the  body  surface. 


4.3.2  Input  Hinge  Points 


The  user  may,  If  the  offset  centerline  geometric  option  Is  not 
used.  Input  the  hinge  point  locations.  This  option  Is  exercised  by  setting 
IHINGE  *  1  and  specifying  values  for  H1(I,L),  the  hinge  point  locations 
In  (x,y)  space  for  each  of  NPHI  planes.  Note  that  HI  Is  a  complex  array, 
where  HI  ■  x+ly.  JA  must  also  be  specified  with  this  option. 

If  an  axlsymmetrlc  mapping  is  desired,  the  parameter  IMAP  Is  set 
to  0  and  only  one  plane  of  HI  data  need  be  Input.  (An  axlsymmetrlc  mapping 
is  assumed  If  the  body  geometry  Is  axlsymmetrlc. )  For  a  non-symmetrl c 
mapping,  NPHI  (see  Section  4.3.3)  planes  of  HI  data  must  be  specified, 
and  each  plane  must  contain  JA+2  hinge  points.  The  HI  points  are  input  in 
the  order  shown  in  Figure  4.7. 

4.3.3  Definition  of  Computational  Grid 

Calculations  in  CM3DT  are  performed  In  a  computational  grid 
constructed  in  the  transformed  (£,n,8)  space,  as  described  in  Reference  1. 
The  coordinate  surfaces  are  then  defined  by  9  a  constant  surfaces  (<f>  planes 
£  *  constant  surfaces  (which  are  nearly  normal  to  the  body),  and  r)  *  con¬ 
stant  surfaces  (which  are  nearly  parallel  to  the  body).  Note  that  these 
surfaces,  when  transformed  back  to  the  physical  (x,y,1>)  space,  result  in 
curved  coordinate  surfaces,  as  illustrated  In  Figure  4.8. 

The  number  of  grid  points  used  in  the  calculation  are  specified 
by  the  parameters  IMAX,  MMAX,  and  NMAX.  LMAX  meridional  planes  are 
equally  spaced  circumferentially  on  the  body,  MMAX  surfaces  are  equally 
spaced  In  £  within  each  meridional  plane  along  the  body  surface,  and  NMAX 
points  are  equally  spaced  In  n  between  the  body  and  shock  along  each 
£  *  constant  curve  within  each  meridional  plane.  The  maximum  allowable 
values  for  LMAX,  MMAX,  and  NMAX  are  9,  18,  and  11,  respectively. 


The  parameter  ISA  Is  used  to  denote  those  cases  In  which  a  pitch 
plane  of  symmetry  exists.  For  ISA  «  1,  the  body  geometry  must  be  symmetric 
about  the  pitch  plane,  there  must  be  no  sideslip  (BETA  •  0.0),  the  angle 
62  must  be  zero,  and,  If  NPHI<0,  there  must  be  no  geometric  offset  of  the 
cross-sections  In  the  yaw  plane  (02  *  0);  all  other  cases  require  ISA  ■  2. 

With  a  pitch  plane  of  symmetry,  calculations  are  performed  only  between 
0  ■  0°  and  <p  *  180*;  with  ISA  ■  2,  $  varies  between  0*  and  360°. 

For  those  cases  with  ISA  -  1,  LMAX  must  be  specified  as  an  odd 
Integer;  with  ISA  ■  2,  LMAX  must  be  an  Integer  multiple  of  4.  (These  re¬ 
strictions  arise  from  the  numerical  treatment  of  the  centerline.)  The  CM3DT 
program  will  automatically  correct  the  value  of  LMAX  specified  to  satisfy 
these  criteria  If  necessary.  However,  If  the  corrected  value  of  LMAX  ex¬ 
ceeds  the  maximum  allowed  (9),  the  program  will  stop. 

4.4  OUTPUT  CONTROLS 

In  the  CM3DT  code,  complete  field  output  Is  provided  every  K0UT 
steps,  starting  with  step  0,  as  well  as  at  the  last  time  step.  The  amount  of 
output  provided  Is  determined  by  the  parameter  JOUT.  For  JOUT  *  0,  only  the 
primary  flow  variables  at  each  grid  point  are  printed  (p,u,s/R,etc.),  as  well 
as  Information  regarding  the  position  and  velocity  of  shock  points  and  the 
shock  slopes. 

When  JOUT  *  1,  the  output  will  also  include  parameters  of  the  co- 
ordinate  transformation  (£,n,g,C,S,$) .  For  non-symmetrlc  mappings  (IMAP  ■  1), 
information  on  the  circumferential  variations  of  the  mapping  at  each  grid 
point  (g^.C^)  will  be  printed  for  JOUT  ■  2.  Note  that  specification  of  an 
ax 1  symmetric  mapping  (IMAP  »  0)  precludes  exercising  the  JOUT  »  2  option, 
since  all  circumferential  derivatives  of  the  transformation  would  then  vanish. 

A  complete  description  of  the  variables  output  using  these  various 
options  may  be  found  In  Section  5.0. 
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In  addition  to  the  printed  output,  CM3DT  can  also  generate  output 
on  a  binary  file.  The  binary  file  output  Is  necessary  for  restarts  (Sec¬ 
tion  4.5.1)  and  for  generating  the  required  Initial  data  for  supersonic 
afterbody  calculations  (Section  6.0).  The  parameter  controlling  the  genera¬ 
tion  of  binary  file  output  Is  KTAPE.  No  binary  file  output  is  created  if 
KTAPE<0;  otherwise,  complete  flow  field  data  Is  written  on  the  file  every 
KTAPE  time  steps. 

4.5  OTHER  PROGRAM  CONTROLS 

4.5.1  Time  Step  Controls 

The  number  of  time  steps  to  De  taken  by  the  CM3DT  code  is  governed 
by  the  inputs  KSTART  and  KMAX.  KSTART  denotes  the  number  of  the  time  step 
at  which  the  calculation  begins;  -If  KSTART  Is  not  input  the  solution  Is 
automatically  Initialized  at  step  0,  and  a  value  of  KSTART  greater  than 
zero  Indicates  a  restart  of  a  previous  solution  from  binary  file  output. 

For  a  restart,  KSTART  must  correspond  to  one  of  the  time  steps  stored  on 
the  binary  file. 

The  last  time  step  to  be  taken  Is  denoted  by  KMAX.  When  this 
step  Is  reached,  the  solution  Is  terminated,  the  field  data  at  that  step 
is  output,  and  (If  KTAPE>0)  the  solution  is  written  on  the  binary  file. 

The  size  of  the  time  step  Is  controlled  internally  by  the  CM3DT 
program.  The  computed  time  step  may  be  altered  by  the  factor  STAB  at  the 
discretion  of  the  user;  however,  values  of  STAB  greater  than  1.0  may  lead 
to  unstable  calculations.  In  normal  applications,  the  user  need  not  be 
concerned  with  modifying  the  default  value  of  this  parameter  (1.0). 

During  the  calculation,  the  computational  scheme  requires  that 
after  each  time  step  an  inverse  mapping  (from  (E,n,9)  space  to  (x,y,4>)space) 
be  performed  to  update  the  transformation  parameters  and  relocate  the  grid 
points  in  physical  space  as  the  grid  varies  with  time.  For  increased 
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efficiency,  however,  this  Inverse  mapping  need  not  be  performed  after  every 
time  step.  The  Input  parameter  KGRID  (preset  to  10)  defines  the  number  of 
time  steps  to  be  taken  between  calls  to  the  Inverse  mapping  subroutine. 

4.5.2  Solution  Initialization 

When  starting  a  solution  from  step  0,  the  CM3DT  code  automatically 
generates  an  estimate  of  the  complete  flow  field  to  provide  the  required 
Initial  data.  Only  one  user  Input  Is  required  for  this  Initialization,  the 
effective  nose  radius,  RN.  The  value  of  RN  controls  the  estimate  cf  the 
shock  stand-off  distance  at  the  centerline,  which  Is  based  on  a  correlation 
of  stand-off  distances  for  spheres  vs.  Mach  number.  Thus,  RN  should  re¬ 
present  the  effective  radius  of  the  nosetlp  being  computed  at  the  centerline. 

4.5.3  A-Scheme  Calculations 

Two  calculation  procedures  are  available  within  the  CM3DT  program: 
the  non-conservation  code  (I LAM  *  0)  and  the  A -differencing  code  ( ILAM  -  1). 
Because  of  its  slightly  greater  efficiency,  the  non-conservation  option  is 
the  preferred  procedure  for  convex  body  shapes  (no  embedded  shocks)  without 
sharp  corners.  For  indented  shapes,  where  embedded  shocks  are  anticipated, 
the  A-differencIng  procedure  Is  required.  Furthermore,  the  A-scheme  offers 
the  advantage  of  minimizing  oscillations  in  the  solution  when  corners  are 
present  in  the  body  geometry  that  are  not  adequately  resolved  with  the  com¬ 
putational  grid. 

4.5.4  Shock  Tables 

In  equilibrium  real  gas  calculations,  downstream  shock  point  pro¬ 
perties  are  determined  from  Interpolation  on  a  table  of  shock  properties  as 
functions  of  the  normal  component  of  freestream  velocity  relative  to  the 
moving  shock.  The  number  of  entries  In  the  table  can  be  specified  through 
the  Input  parameter  NUWTB,  which  is  preset  to  91.  The  upper  and  lower  limits 
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on  the  table  are  specified  as  UWTBMX  and  UWTBMN,  respectively,  and  are 
preset  to  1.2  and  0.1.  (These  values  represent  the  extremes  In  the  ratio 
of  the  normal  velocity  component  relative  to  the  shock  to  the  freestream 
velocity. ) 
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SECTION  5 
CODE  OUTPUT 


Standard  output  of  the  CM3DT  code  includes  printout  of  the 
Input  variables  and  body  geometry  data,  parameters  of  the  coordinate 
transformation,  a  brief  time  history  of  the  calculation,  and  complete 
field  data  every  KOUT  steps.  Optional  outputs  are  available  If  requested 
by  the  user,  as  described  below.  A  glossary  of  output  variables  is 
provided  in  the  Ippendix  (Section  10. 2 ), and  examples  of  the  code  output 
are  also  provided  in  the  Appendix  (Section  10.5) . 

5.1  PRELIMINARY  OUTPUT 

5.1.1  Input  Variables 

The  first  page  of  output  generated  by  the  CM3DT  code  lists 
values  of  input  variables  defined  by  the  user,  as  well  as  information 
computed  Internally  by  the  code  based  on  the  Input  parameters.  First, 
information  on  all  of  the  freestream  properties  Is  provided,  followed  by 
data  concerning  the  thermodynamic  model  (ideal  or  real  gas)  selected  for 
the  calculation. 

Information  is  then  provided  about  the  angle  of  attack  and 
yaw  angle  relative  to  the  nosetip  centerline  as  computed  by  the  code 
from  the  input  angles  of  attack  and  sideslip  (which  are  defined  relative 
to  the  afterbody  axis).  Following  this,  the  input  variables  defining  the 
computational  coordinate  arid,  the  number  of  time  steps  to  be  taken,  and 
the  output  controls  are  listed. 


This  first  page  of  output  also  provides  values  of  the  theoretical 
stagnation  pressure  (PSTAG)  and  total  enthalpy  (HTOT),  as  well  as  values  of 
normal  shock  pressure  (PSH),  entropy  (SRSH),  and  Mach  number  (AMSH).  as 
computed  Internally  by  the  code  from  the  specified  freestream  conditions. 

Finally,  this  page  of  output  concludes  with  explicit  statements 
of  the  options  that  have  been  specified  by  the  user, 

5.1.2  Body  Geometry  Definition 

The  next  section  of  the  output  details  the  specification  of  the 
body  geometry  for  the  calculation.  For  an  analytic,  axl symmetric  geometry, 
the  values  of  the  coefficients  (ACOF)  of  each  curve  segment  are  printed* 
as  well  as  the  x-values  of  each  segment  end  point. 

For  arbitrary  body  geometries  without  an  offset  centerline,  the 
tabulated  values  of  XB  and  YB  In  each  PHI  plane  specified  by  the  user  are 
printed,  followed  by  a  statement  of  which  arbitrary  body  geometry  option, 
as  described  in  Section  4.1.2,  Is  used. 

For  arbitrary  body  geometries  with  an  offset  centerline,  the  values 
of  the  arrays  Dl,  D2,  AE,  BE,  CE,  and  DE  are  printed  out  at  each  axial  station, 
XE(I).  Information  Is  also  printed  out  defining  the  orientation  of  the  new 
computational  nosetlp  axis  relative  to  the  afterbody  axis,  as  well  as  the 
values  of  a  and  B  relative  to  the  computational  axis,  as  computed  auto¬ 
matically  within  CM3DT. 

5.1.3  Conformal  Transformations 

Following  the  body  geometry  information,  the  output  provides  in¬ 
formation  on  the  hinge  points  (whether  generated  automatically  by  the  code 
or  Input  by  the  user)  In  the  physical  (x,y,<f>)  space,  as  well  as  the  hinge  • 
point  locations  after  each  Intermediate  conformal  transformation.  The  ex¬ 
ponents  (POWER)  of  each  Intermediate  transformation  in  each  meridional  plane 
are  also  provided.  Note  that  the  hinge  points  are  printed  as  H(I,J,L),  a 
complex  number,  where  the  subscripts  identify  the  I—  hinge  point  after 
the  (J-l)—  transformation  in  the  L~  meridional  plane. 
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The  final  portion  of  this  section  of  output  defines  the 
stretching  parameters  (AL)  used  In  the  final  transformation  In  each 
meridional  plane,  as  well  as  Information  defining  the  downstream  boundary 
of  the  region  being  computed  In  the  transformed  space  (XIL,  XILTH) . 

5.1.4  Transformed  Body  Geometry  Information 

The  last  part  of  the  preliminary  output  of  the  CM3DT  code 
provides  body  geometry  Information  at  each  computational  grid  point 
located  on  the  body  surface.  This  Information  Is  output  for  each  meridional 
plane  being  computed,  and  Includes  data  on  the  location  of  the  point  and 
the  geometric  derivatives  In  both  the  physical  space  (XB,  YB,  DYDX,  DYDPH) 
and  in  the  transformed  space  (XIB,  ETAB,  ETABXI,  ETABTH). 

5.2  STANDARD  OUTPUT 

5.2.1  Field  Data 

Complete  field  data  Is  output  by  the  CM3DT  code  every  KOUT 
time  steps,  starting  with  step  0,  as  well  as  at  the  final  time  step,  KMAX. 
Within  every  meridional  plane,  the  data  output  consists  of  the  location 
of  every  grid  point  In  physical  space  (X,Y),  the  local  static  pressure  (P), 
entropy  (SR),  density  (RHO),  and  Mach  number  (MACH),  and  the  velocity 
components  In  (£,n,0) space  (U,V,W).  In  addition,  the  variation  between 
the  computed  and  theoretical  values  of  total  enthalpy,  ratioed  to  the 
known  theoretical  total  enthalpy,  is  provided  as  a  measure  of  convergence 
(DH). 


Information  is  also  provided  within  each  meridional  plane  on 
the  shock  position,  velocity,  and  slopes  In  the  transformed  space.  The 
shock  position  is  defined  by  XIS  and  ETAS,  and  the  shock  velocity  In  the 
transformed  space  along £»  constant,  6  «  constant  curves  is  given  as  ETAST. 
This  transformed  shock  velocity  corresponds  to  physical  shock  velocity 
ETSTOG.  Shock  slopes  in  the  transformed  space  are  given  by  ETASXI  and  ETASTH. 
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5.2.2  Time  History  Output 


To  provide  a  history  of  the  time-dependent  calculation,  the 
CM3DT  code  prints  at  every  time  step  values  of  four  quantities  that 
provide  a  means  of  judging  the  convergence  of  the  solution.  (Convergence 
criteria  for  this  code  are  discussed  In  Section  6  of  Volume  I  of  this 
report.)  The  four  quantities  printed  are  PSTA6  (pressure  at  the  center¬ 
line  body  point).  DELTA  (distance  between  the  body  and  the  bow  shock  at 
the  centerline),  WRMS  (root-mean-square  of  the  shock  velocities),  and 
ETSOGM  (the  maximum  shock  velocity).  The  convergence  of  PSTA6  and  DELTA 
Indicate  the  potential  convergence  of  the  calculation  In  the  vicinity  of 
the  centerline,  while  WRMS  and  ESTOGM  serve  as  Indications  of  the  asymptotic 
approach  of  the  entire  solution  to  the  steady  state  limit. 

5.3  OPTIONAL  OUTPUT 

5.3.1  Optional  Field  Output 

Optional  field  output  will  be  generated  by  the  CM3DT  code 
when  the  Input  parameter  JOUT  Is  set  to  1  or  2.  This  optional  field  out¬ 
put  Is  printed  Immediately  following  the  standard  field  output,  every 
KOUT  time  steps. 

With  JOUT  »  1,  information  Is  provided  on  the  conformal 
transformations  within  each  meridional  plane.  Specifically,  at  every  grid 
point,  values  are  printed  for  the  (E,n)  location  of  the  point  in  trans¬ 
formed  space  (XI, ETA),  and  Information  on  derivatives  of  the  transforma¬ 
tion  within  that  meridional  plane.  These  derivatives  are  g  ■  ■  Ge^t0 

and  *  *  3(1  og  g)/3c,  and  outputs  are  provided  for  the  magnitude^  g  (AG), 
the  real  and  imaginary  parts  of  g  (REG,  IMG),  the  real  and  imaginary  parts 
of  e’iu(CISl,  CIS2),  and  the  real  and  imaginary  parts  of  $  (PHI1,  PHI2). 
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If  JOUT  *  2,  the  same  Information  provided  with  OOUT  *11$ 
printed,  as  well  as  Information  concerning  the  circumferential  derivatives 
of  the  transformation.  Values  are  printed  for  the  reel  and  Imaginary 
parts  of  d;/d$  (REZETPHI,  IMZETPHI)  and  for  the  real  and  Imaginary  parts 
of  3g/d$  (REGPHI ,  IMGPHI)  at  every  computational  grid  point  when  this 
option  Is  exercised.  Note  that  the  JOUT  ■  2  option  cannot  be  Invoked  If 
an  axlsyimetrlc  coordinate  mapping  (IMAP  ■  0)  has  been  specified,  since 
these  derivatives  then  vanish. 

5.3.2  Shock  Table  Output 

If  a  real  gas  calculation  is  being  performed  (IRG  ■  1),  the 
user  may  have  the  shock  property  tables  printed  by  specifying  ITBOUT  *  1. 
When  this  option  Is  selected,  values  of  pressure  (P),  entropy  (S/R),  the 
1  sen tropic  exponent  (G),  and  the  downstream  normal  component  of  velocity 
relative  to  the  shock  (UW2)  are  printed  out  as  functions  of  the  upstream 
normal  component  of  velocity  relative  to  the  shock  (UW1). 
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SECTION  6 


COUPLING  OF  NOSETIP  CODE  TO  SUPERSONIC  AFTERBODY  CODE 


The  CM3DT  nosetip  flow  field  code  has  been  coupled  to  the 
supersonic  afterbody  code  described  in  References  2  and  3  to  provide  a 
capability  for  the  prediction  of  total  vehicle  aerodynamics  for  reentry 
vehicles  with  ablated  nosetips.  This  coupling  requires  interpolation  on 
the  nosetip  solution  to  provide  the  required  initial  flow  information  at 
the  starting  line  of  the  afterbody  code.  The  relationship  between  the  nose¬ 
tip  computational  (x ,y ,d>)  and  afterbody  (r,e,z)  cylindrical  coordinate 
systems  is  described  below,  as  well  as  the  inputs  required  to  couple  CM3DT 
to  the  afterbody  solution. 

As  discussed  in  Section  4,  the  nosetip  computational  axis  may  be 
rotated  from  the  afterbody  axis  in  both  pitch  and  yaw,  the  angles  of  rota¬ 
tion  being  denoted  by  ^  (DELI)  and  5^  (DEL2),  respectively.  (Note  that  if 
the  computational  axis  differs  from  the  geometric  definition  axis,  the  values 
of  and  62  are  recomputed  within  CM3DT,  as  described  in  the  Appendix, 
Section  10.4.)  The  angles  of  attack  (a)  and  sideslip  (B)  required  as  input 
to  the  CM3DT  code  are  defined  relative  to  the  unbent  afterbody  axis. 

For  the  nosetip  calculation,  however,  CM3DT  requires  the 
effective  angles  of  attack  (0^)  and  sideslip  (3^)  relative  to  the  nosetip 
centerline.  These  quantities  are  computed  internally  by  CM3DT  from 


sin  8n  =  sin  B  cos  <$2  +  cos  a  cos  6  s ™  $2  cos  5A 
-  sin  a  cos  B  sin  S2  sin  6^ 


(6.1) 


tan  aN  =  Vy/Vx 


(6.2) 


..wafer:  .vuMfept”*' 


where 


vy  *  Vco  tcos  Ot  cos  6  (sin  61  cos  $2  cos  6A  +  cos  $1  sin  6ft) 

+  sin  a  cos  3  (-sin  6j  cos  6g  sin  +  cos  cos  $A) 

-  sin  6  sin  sin  621  (6.3) 

Vx  s  Vro  [cos  a  cos  3  (cos  6j  cos  6g  cos  6A  -  sin  sin  6A) 

+  sin  a  cos  B  (-cos  6^^  cos  sin  6A  -  sin  61  cos  6^) 

-  sin  8  cos  sin  fig]  .  (6.4) 

It  is  required  that  the  nosetip  computational  axis  intersect  the 
afterbody  axis;  this  point  of  intersection  has  an  x  value  of  xQ  (X0)*1n 
the  nosetip  (x,y,<j>)  coordinate  system  and  a  z  value  of  zQ  (ZO)  in  the 
afterbody  (r,e,z)  coordinate  system.  If  the  afterbody  axis  is  bent,  the 
bend  occurs  at  a  z  value  of  z  (ZB).  These  parameters  are  illustrated  in 
Figure  4.1,  and  must  be  included  in  the  input  to  the  afterbody  code. 

The  relationship  between  the  afterbody  (r,0,z)  cylindrical  co- 
and  the  (x,y,4>)  nosetip  computational  cylindrical  coordinate  system  Is 
described  in  the  Appendix,  Section  10.4. 

In  the  generation  of  starting  line  data  for  the  afterbody  code, 
it  Is  necessary  to  completely  define  the  shock  layer  flow  in  a  z  =  constant 
surface,  at  grid  points  located  at  known  0  values  at  known  percentages  of 
the  shock  layer  thickness.  The  first  step  in  generating  the  starting  line 
data  is  to  locate  the  shock  position  in  the  initial  data  plane;  i.e., 
r  =  rg  (0,z)  for  the  known  values  of  z  and  0.  This  requires  an  iterative 


♦The  parameter  XO  In  the  afterbody  code  is  identical  to  the  parameter  XHO 
in  the  CM3DT  code. 
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procedure,  assuming  a  value  of  r$,  determining  the  corresponding  values  of 
and  checking  the  value  of  y  so  determined  against  the  shock  position 
m  the  nosetip  coordinate  system,  y  «  y.  (Xi*).  Once  the  shock  position  In 
the  Initial  data  plane  has  been  determined,  the  location  of  all  afterbody 
grid  points  in  the  Initial  data  plane  Is  known. 


Knowing  (r.e.z),  the  corresponding  value  of  (x,y,<p)  can  be  cal¬ 
culated,  and  the  flow  variables  P,s,u',v'f  and  w'  can  be  determined  by 
interpolation  on  the  nosetip  solution.  The  nosetip  transformed  velocity 
components  (u'.v’.w')  may  be  used  to  compute  the  nosetip  cylindrical 
velocity  components  (U,V,W)  from 

U  =  u'C  -  v'S 
V  ■  u'S  +  v'C 
W  *  W 


(6.5) 

(6.6) 


(where  C  +  IS  *  e  iw).  The  afterbody  cylindrical  vel 
u  (radial)  and  v  (circumferential)  are  computed  from 


ocity  components 


u  =  cos  e  vy  +  sin  e  V2 
v  *  -sin  9  V  +  cos  9  V, 

•J  Z 


(6.8) 

(6.9) 


where 


vy  35  U(-sm  6a  cos  d1  cos  6Z  -  cos  6A  sin  6j) 

+  ^  C0S *  ~  W  sir,*H-s1.n  6a  sin  cos  6^  +  cos  $A  cos  6j) 
-  (V  s1n4>  +  W  cos<1>)  sin  6A  sin  $2 

V2  *  “’Ll  cos  61  sin  62  -  (v  cos*  -  W  s1nd»)  sin  6j  sin  «2 
+  (V  si nd>  +  W  cos*)  cos  62 


(6.11) 


The  afterbody  axial  component  w  can  be  determined  from  the  conservation  of 
total  enthalpy,  expressed  as 

H  «  h(p,s)  +  ^  (u2  +  v2  +  w2)  (6.12) 

where  H  is  the  total  enthalpy,  which  is  constant  and  equal  to  the  free- 
stream  total  enthalpy  for  a  steady  inviscid  flow. 

Forces  and  moments  on  the  nosetip  are  determined  by  integration 
in  the  nosetip  coordinate  system,  referencing  the  moments  to  the  center  of 
gravity  location  specified  In  the  afterbody  code  by  (X,Y,Z),  as  shown  in 
Figure  6.1.  The  relationship  between  the  nosetip  computational  and  after¬ 
body  coordinate  systems  is  detailed  in  Section  10.4*,  (iu.yg.Zg)  represents 
the  center  of  gravity  location  In  the  nosetip  coordinate  system. 

The  force  integrals,  evaluated  relative  to  the  nosetip  computational 
axes  and  using  the  sign  conventions  depicted  in  Figure  6.2,  may  be  expressed 
as 


2ir  f  Xjr(<i>) 

pVb) 

*B0 


n 


dxdcp 


xe(4>) 

pyb(cos<j>  +  sin4>  yb$/yb)dxd<t> 
o 


xE(d>) 

pyb(sin<J>  -  cost})  yb<(/yb)dxd<i> 
o 


(6.13) 


(6.14) 


(6.15) 
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Figure  6.2,  Force  and  Moment  Sign  Conventions 


where  Fan,  Fnn,  and  Fyf)  represent  the  axial,  normal,  and  side  forces  on 
the  nosetlp,  respectively.  In  these  expressions,  p  Is  the  surface 
pressure,  xb0  Is  the  x  value  of  the  centerline  body  point,  and  Xg(d> )  Is 
the  x  value  at  the  body  surface  In  the  afterbody  Initial  data  plane, 
expressed  In  the  nosetlp  coordinate  system.  The  nosetlp  moment  Integrals 
may  be  expressed  as 


ft 

O  Xg 


xE(<J>) 

pybtx(cos<|>  +  sin*  y^/y 
o  +  ybybx  cos<pi  d*d<t> 


b> 


+  <xo  +  x2)FNn  -  *2FA 


N 


xE(<j>) 

pyb[x(sin<j>  -  cos *  yb(j)/yb 

0  +  Vbx  s1n<l>1  dxd<t> 

+  (x  +  x«)Fy  +  ip  fa 

0  z  Yn  i  An 


(6.16) 


(6.17) 


xE(<J>) 

pybz[2sin(f)  coscp  +  ( si n2(p  -  cos 2<J> ) 
0  VV  dxd»  -  J2FNn  • 


(6.18) 


where  M^,  MyN,  and  MyN  represent  the  nosetlp  pitching,  yawing,  and 
rolling  momrnts,  respectively,  using  the  sign  conventions  depicted  in 
Figure  6.2,  with  x2,  y2,  and  z 2  defining  the  center  of  gravity  location 
in  the  nosetlp  Cartesian  coordinate  system. 
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These  nosetlp  forces  end  moments  cen  be  expressed 
to  the  afterbody  coordinate  system  as 

Fa'  «  (cos  cos  cos  6a  -  sin  sin  5A)FA 

+  (sin  cos  <52  cos  6^ 

+  cos  sin  5A)FN  -  sin  62  cos  5aFYn 


*  -  (cos  cos  sin  5a  +  s^n  61  cos  ‘V^A 

-  (sin  6,  cos  $2  s-in  5a 

-  cos  5j  cos  6A)FN  +  sin  ^ ' 


A  Y. 


Fu  =  cos  5,  sin  SoF,  +  sin  6,  sin  60F, 


Vm  2  A, 


N 


1  2  Nk 


+  cos  ^Fy 


M-7 1  =  -cos  5,  sin  <5«MV  -  sin  5,  sin  5,MV 
Z  1  nN  l  i  rN 

+  cos  59M7 
2  2N 


M 


Y 


I 


=  -  (cos  <5j  cos  sin  <5A  +  sin  6 

-  (sin  6^  cos  sin  6A  -  cos  5^ 

-  sin  5,  sin  <5flM7 

t  A 


j  cos  6a)M 

cos  6. )MV 
a  tn 


MX 


*  (cos  6^  cos  cos  6a  -  sin  6j  sin  $A)MX 

+  (sin  5^  cos  6g  cos  5a  +  cos  51  s^n  5A^MY 

+  sin  6,  cos  5«M7 
Z  a  zN 


N 

N 


relative 


(6.19) 
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(6.21) 


(6.22) 


(6.23) 


(6.24) 
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The  axial  force  (F^')»  pitching  moment  ( ' ) .  and  yawing  moment  (My')  of 
the  cosetip  in  the  afterbody  coordinates  are  denoted  by  primes  to  Indicate 
that  these  quantities  have  been  computed  assuming  zero  base  pressure 
(pfa  *  0).  In  the  afterbody  code  to  which  CM3DT  has  been  coupled,  however, 
forces  and  moments  are  computed  assuming  that  pb  *  p^.  For 
consistency,  the  axial  force  and  yawing  and  pitching  moments  for  the 
nosetip,  expressed  in  the  afterbody  coordinates,  can  be  modified  to 
fit  this  assumption  using 


fa 

s  V 

•ip, 

L 

2tt 

/  rb*(e)de 

(6.25) 

V 

Mz 

=  Y 

’  Poo 

*0 

rzv 

/  i^rbcos9  -  |(Y-d)]rb2(e)de 

(6.26) 

► 

My 

*  My' 

+  p. 

J 

0 

/■2tt 

I  <!rb sine  +  1  ^)rb2(0)de 

(6.27) 

J 

0 

where  the  integrals  are  to  be  evaluated  at  the  initial  data  plane  of  the 
afterbody,  in  the  afterbody  cylindrical  coordinate  system  (r,9,z). 


i 

i 
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SECTION  7 


INPUT  GUIDE  LIST  AND  SAMPLE  CASES 


7.1  INPUT  GUIDE  (CM3DT) 

The  first  input  record  is  a  title  card  which  may  contain  up  to 
80  characters  of  alphanumeric  information.  Following  the  title  card, 
the  input  variables  are  read  in  the  NAMELIST  format,  where  the  NAMELIST 
group  name  is  INPUT. 


INPUT  CHECKLIST 

(Title  Card;  Alphanumeric  information 
in  columns  1-80) 

SINPUT  (Namelist  Identifier) 

(Freestream  Conditions) 

AMINF 

QINF 

IRG 

IATMP 

ALT 

PINF 

TINF 

G I  DEAL  \ 

RIDEAL 
Z I  DEAL 
ALPHA 
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BETA 


DELA 

DELI 

DEL2 

NS  EG 

ACOF(I.J) 

ADIV(J) 

NPHI 

PHIBD(L) 

NPTS(L) 

XB(I.L) 

YB(I.L) 

NPHI 
XE  ( I ) 

AE(  I ) 

BE  ( I ) 

CE  ( I ) 

OE(I) 

D1(I) 

D2(I) 

NPTS(l) 

NSEG 


(Geometry  -  Analytic) 


(Geometry  -  Arbitrary,  without  Offset) 


(Geometry  -  Arbitrary,  with  Offset) 


/ 
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ACOF(I,J) 

AOIV(J)  ■  4 

XHO 

(Grid  Definition) 
XBO 

XBB(L) 

NMAX 

MMAX 

LMAX 

IMAP 

IHINGE 

Hl(I.L) 

DELH 

JA 

ISA 

(Time-Step  Controls) 

KMAX 

STAB 

KSTART 

KGRID 

(Output  Controls) 

KOUT 

JOUT 

KTAPE 

ITBOUT 
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(Miscellaneous  Controls) 

I  LAM 
RN 

TEST ( I ) 

(Shock  Table  Controls) 
(Real  gas  only) 

NUWTB 

UWTBMX 

UWTBMN 

(End  of  NAMELIST  Input) 


7.2  INPUT  GUIDE  (AFTERBODY  CODE) 


The  standard  inputs  to  the  afterbody  code  are  defined  in  Reference 
4.  The  version  of  this  afterbody  code  modified  for  use  with  CM3DT  has  four 
additional  input  parameters  that  must  be  defined  to  relate  the  nosetip 
and  afterbody  coordinate  systems.  These  inputs  are: 


ZB  z.  Afterbody  axial  location  where 

afterbody  axis  is  bent 

ZO  z0,  Afterbody  axial  location  where 

afterbody  and  nosetip  axes  intersect 

XO  x0,  Nosetip  computational  coordinate  system 

axial  location  where  afterbody  and  nose¬ 
tip  axes  intersect  (not  required  if  off¬ 
set  centerline  nosetip  geometry  option  used) 

TEST (4)  Convergence  criterion  for  shock  position 

interpolation  in  generating  starting  line 
data,  preset  to  0.001. 


7.3 


SAMPLE  CASES 


In  this  section  the  Inputs  required  for  three  sample  calcula¬ 
tions  with  the  CM3DT  code  are  presented,  along  with  the  Inputs  required 
for  a  sample  afterbody  calculation.  The  samples  of  output  presented  In 
the  Appendix,  Section  10.5,  have  been  selected  from  these  sample  cases. 


Case  1  -  The  first  sample  case  Is  a  sphere  In  Ideal  gas  (y  *  1.4) 
at  M»  ■  20  in  both  pitch  and  yaw,  with  a  ■  5°  and  S  *  5°.  The  geometry 
Is  defined  analytically,  and  an  axlsynmetrlc  coordinate  mapping  Is  used. 

The  non-conservation  version  of  the  code  Is  used  for  this  calculation. 


INPUT  FOR  CASE  1 


CM3DT <NC>  SPHERE  MftCH  £0  PLPHft*5  BETft=3 
$ INPUT 

AMINF*£0.  0»  ALPHR*5.  0i»  BETfi*5.  0» 

NSES»1  f  HCOF l)«i.  0»0.5'  0.5* 

KB 0*0.  0»  XBB  <’l  >  *0.  5»  IMPP*0  j  JR*5 1> 

MMPX«1 0>LnAX«£» I SP*£> 

KMflX«300» KQUT«300» 

ILftM*0»  RN*0. 5»  3 


Case  2  -  The  second  case  Is  a  calculation  of  an  Ideal  gas  flow 
at  *  11.6  over  an  axlsymmetric  Indented  shape,  the  Mildly  Indented  Body 
(MIB),  depicted  in  Figure  7.1  and  described  in  Reference  5,  at  a  *  2°.  An 
arbitrary,  axisymmetrlc  geometry  definition  is  used,  and  an  axlsymmetric 
coordinate  transformation  is  Invoked.  Since  the  MIB  shape  is  Indented, 
with  an  embedded  shock,  the  X-dlfferencing  calculation  procedure  is  used. 
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INPUT  FOR  CASE  2 


CM3DT  (LAMBDA)  MIB  MACH  11.6  ALPHA=2 
SINPUT 

AMINF«11.6?ALPHA»2. 0? 

MPHI *1  ? NPTS  <1>  *30? 

XB<1>  *0.  ? .  0065? .  0269!. .  0617? .  1096?.  1701  ? .  2672? .  4213?  .  61 1 3? .  3154? 
1 . 0332?  1. 291 2 ? 1.5692*  1 .3542*2. 1395?  2. 4122;.  2.  6584?  2.  8943? 3.  1255? 
3. 3684?  3. 6552?  4.  0?4.2?4.4?4.6?4.3?5.0?5.5?6.0?6.5? 

VB  <1 > *0. ? . 1511 ? .  3059* . 4606? . 6097? . 751 ? . 379? . 9803?  1. 07 27? 1 . 1722? 
1 .2809 ? 1.4043!*  1.551 1»  1. 7299?  1. 942?2. 1818?  2.  431 3. •  2. 6766?  2.  £076? 
3.  0934s.  3.2645?  3.  3935!.  3.  4433*  3.  4793s.  3.  5004s*  3. 5214;.  3.  5424?  3.  595? 
3, 6475?  3. 7001 ? 

XB0=O.  0?  XBB  <1>  *9*6.  0»  IMAP=0? JA=7? 

NMRX=11?MMAX*18?LMAX*9? 

K  MAX  * 5  0  0  ?  K  OU  T  *5  0  0  ? 

I LRM=1 ?  RN=1 .  75?  $ 


Case  3  -  The  last  sample  case  for  the  CM3DT  code  has  been  selected 
to  demonstrate  the  use  of  an  asymmetric  conformal  mapping  procedure.  The 
nosetip  being  computed  is  the  Blunt-1  asymmetric  shape,  with  offset  circular 
cross-sections  assumed  (NPHI=2).  The  calculation  is  performed  with  the- 
A-differencing  procedure,  at  a  Mach  number  of  13.4  in  ideal  gas,  at  a  »3°. 
Note  that  the  nosetip  centerline  is  rotated  in  the  pitch  plane  by  5^=  9.648° 
from  the  afterbody  axis,  leading  to  an  effective  nosetip  angle  of  attack  of 
aN  =  12.648°.  This  nosetip  shape,  and  Its  orientation  relative  to  the  after¬ 
body  axis,  is  illustrated  in  Figure  7.2. 
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INPUT  FOR  CASE  3 

CM3DT (LAMBDA)  BUJNT-1  NOSE  SHAPE  MACH  13.4  AlPHA=3 

* INPUT 

AMINF=l 3. 4 ? ALPHA*3. Ok  DEL  1=9. 648 ? 

NPH I =2 ?  NPT 3  < 1 > =28 ?  3 0 ? 

XB  t.l  >  \  )  22? .  23?  .247? .  28?  .3?  .34?  -4s.  .427? .  5? .  7? .  3? .  9?  1 . 0?  1 .  1  ?  1 .  IS?  1 . 2 
1 • 22? 1 . 24? 1 . 28? 1.4? 1.6? 1 . 8?  2. 0?  2. 2?  2. 4?  2. 6?  2. 3? 3. 0? 

VB  <1  ?  1 )  =  0.  ? .  2? .  3? .  4? .  433?  .5?.  570?.  6?.  667?.  81?.  365?  .  92?  .  96?  1 . 0?  1 .  022? 
1 . 038?  1 .  043? 1 . 05  ? 1 . 06?  1 . 095? 1 . 1 5?  1 . 203? 1 . 265  ? 1.321?  1.377? 1.433? 

1.489? 1.545? 

XB  v 1 ?  2)  =  .22? . 23  ? . 245  ? . 26  ? . 3  ? . 34  ? . 33  ? . 42  ? . 46  ? . 6  ? . 8  ? 1 . 0  ? 1 . 2  ? 1 . 4  ? 1 . 44  ? 
1.43?  1.53?  1.56'?  1.6?  1.65?  1.7?  1.74?1.S?  1 . 9  ?  2 .  0?2»2?2.4?2.6?2.3?3.0? 

VB  <.  1  ? 2)  =  0.  ? .  04? .  1  ? .  16? .  245? .  302? .  345?  .  378? .  412? .  527? .  695?  .36?  1 . 027? 
1.195? 1.23? 1.26? 1.3? 1. 32? 1.342? 1.363? 1.375? 1.33? 1 . 377? 1 . 372? 1 . 365? 
1.352? 1.339? 1.327? 1.314? 1.301? 

XB0=0. 22?  XBB  <1>  =2. 4?  2. 8? IMAP=1 ?  JA=7? DELH=. 1 5 ? 

NMAX*6?  MMAX=13? LMAX=9? 

KMAX*500?  KOUT=500?  KTAPE-500? J0UT=2? 

ILAM=1?RN=1. 0?* 
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Figure  7.2,  BLUNT-1  Configuration 


Case  4  -  This  sample  case  presents  the  inputs  required  to  run  the 
afterbody  code,  using  the  CM3DT  code  to  generate  the  necessary  Initial  data. 
This  calculation  is  for  a  6°  cone  (with  no  bend  In  the  afterbody  axis), 
with  the  Blunt-1  asymmetric  nosetlp  attached.  The  afterbody  calculation 
Is  performed  at  a  *  3°,  and  the  Initial  data  is  obtained  from  the  CM3DT 
calculation  described  In  Sample  Case  3. 

INPUT  FOR  CASE  4 

1N3IS  BLUNT-1  NOSE  SHAPE  6  DEG  CONE  MACH  13.4  ALPHAS 

ISTAPT»500> 
nmax*h,  ririAX*9» 

MPRINT=0« 

K  0=3 .  4  *  Z  0= 1 3 . 35  ? !» 

ZSTART«12.?5?!.ZEND=4F.  0« 

ALPHA*3!. 

VTAB  a , 1 ) *- 1 , 34 08 » Y TAB  Cl j  £> »-4 . *4 , 

2XYTAB»1£.75?,4?.0, 


SECTION  8 


DEFINITION  OF  APPLICABILITY  OF 
CM3DT  NOSET  IP  CODE 


The  CM3DT  invlscld  nosetlp  flow  field  code  has  extended  the  range 
of  nosetlp  geometries  for  which  successful  Invlscld  flow  field  calculations 
can  be  performed  by  using  a  new  coordinate  system  which  can  be  closely 
aligned  with  the  body  geometry,  and  by  using  the  X-dlfferencIng  scheme  to 
compute  embedded  shocks  on  Indented  shapes  In  an  approximate  manner. 
Limitations  do  exist,  however.  In  the  application  of  this  technique  to 
ablated  reentry  vehicle  nosetips. 

As  an  inviscid  technique,  CM3DT  is  not  capable  of  treating  flows 
where  viscous  effects  are  Important,  such  as  separated  or  oscillating  flows, 
as  might  occur  on  severely  Indented  shapes.  In  addition,  since  CM3DT 
approximates  embedded  shocks  as  isentropic  compressive  discontinuities, 
the  accuracy  of  the  CM3DT  results  diminishes  as  the  strength  of  an  embedded 
shock  grows. 

This  section  details  the  results  of  an  effort  to  define  those 
indented  nosetip  shapes  for  which  the  CM3DT  code  can  be  expected  to  provide 
reliable  flow  field  results.  This  effort  is  restricted  to  indented  nose- 
tip  geometries,  since  CM3DT,  by  virtue  of  its  unique  body-oriented  coordi¬ 
nate  system,  does  not  suffer  from  any  limitations  on  its  application  to 
convex  geometries,  other  than  the  requirement  for  having  a  sufficient 
number  of  mesh  points  to  adequately  resolve  critical  features  of  the  flow 
field.  (This  restriction,  of  course,  applies  to  al_]_  flow  field  techniques.) 

The  approach  taken  to  define  the  applicability  of  the  CM3DT 
technique  to  Indented  nosetip  shapes  was  to  perform  calculations  with  an 
axi symmetric  version  of  the  CM3DT  code  on  parametrically  varying  indented 
nosetip  geometries.  The  generalized  nosetip  geometry  used  in  these 
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calculations,  taken  from  Reference  6,  Is  shown  In  Figure  8.1,  and  Is  defined 
In  terms  of  the  nose  radius  (R^ ) •  the  fillet  radius  (Rp),  the  shoulder  radius 
(R$),  the  cone  angles  ( ©i ,  62).  the  nosetlp  length  (£),  and  the  fineness  ratio 
U/d).  Most  indented  nosetlp  geometries  can  be  modeled  (albeit  approximately) 
with  this  parametric  representation. 

All  calculations  were  performed  with  constant  nose  radius 
(RN  *  0.5);  because  of  the  lack  of  an  Intrinsic  length  scale  in  Invlscid 
flows,  the  effects  of  variation  in  R^  can  be  investigated  by  keeping  Rr 
constant  and  varying  other  parameters.  The  shoulder  radius  was  also  held 
constant  for  these  calculations  (R$  =  0.5),  since  this  parameter  does  not 
significantly  affect  the  flow  in  the  indented  region. 

The  bulk  of  the  calculations  were  performed  with  zero  fillet  radius 
(Rp  a  0),  corresponding  to  a  sharp  compression  corner  in  the  body  geometry. 

A  limited  njmber  of  calculations  were  performed  with  a  non-zero  fillet 
radius,  providing  a  direct  comparison  of  the  effects  of  a  gradual  versus  a 
sharp  compression  on  the  flow  field  computation. 

Most  calculations  were  made  at  H*  =  10,  assuming  ideal  gas  thermo¬ 
dynamics.  Several  calculations  were  made  at  M*  *  5  and  M*.  *  20  in  order  to 
assess  the  effects  of  freestream  Mach  number  on  the  flow  field  results. 

Since  only  indented  shapes  were  of  interest  in  this  study,  the 
configurations  considered  were  constrained  by  the  requirement  that  02>0i. 
Furthermore,  since  embedded  shocks  can  occur  only  when  the  shock  layer  flow 
upstream  of  the  shock  is  supersonic,  values  of  0^  greater  than  40°  were 
not  considered. 

Use  was  made  of  an  axi symmetric  version  of  CM3DT  in  order  to 
minimize  the  computer  resources  required  for  this  effort.  However,  all 
axisymmetric  calculations  were  performed  using  no  more  mesh  points  along 
the  body  (MMAX)  and  between  the  body  and  the  shock  (NMAX)  than  are 
available  within  the  CM3DT  code,  so  that  the  results  obtained  can  be  applied 


57. 


directly  to  the  CM3DT  code.  The  effects  of  angle  of  attack  on  the  appli¬ 
cability  of  the  CM3DT  code  were  assessed  by  examination  of  three-dimen¬ 
sional  calculations  performed  as  part  of  the  validation  procedure  for  the 
CM3DT  technique. 

For  each  calculation,  an  assessment  of  the  quality  of  the 
solution  was  made  in  order  to  determine  if  the  solution  was  acceptable 
(i.e.,  was  the  numerical  technique  capable  of  providing  an  acceptable 
solution  for  that  particular  nosetlp  geometry).  Although  the  evaluation 
of  solution  quality  is  subjective,  several  key  parameters  were  carefully 
examined.  In  particular,  values  of  OH  (variation  of  computea  total  en¬ 
thalpy  from  the  known  theoretical  total  enthalpy)  and.  WRMS  (root-mean- 
square  of  shock  point  velocities)  were  considered;  ideally,  these  quantities 
should  vanish  in  the  limit  of  a  perfectly  converged  solution. 

As  discussed  In  Volume  I  of  this  report  (Reference  1),  criteria 
on  DH  and  WRMS  have  been  developed  to  define  a  converged  solution  for  the 
CM3DT  code.  ’These  criteria  represent  sufficient,  but  npt  necessary, 
conditions  for  convergence.)  In  the  judgement  of  the  suitability  of  the 
solutions  obtained,  these  criteria  were  not  adhered  to  rigidly.  For 
example,  the  formal  convergence  criterion  on  DH  requires  that  ( DH }  not 
exceed  C.05  at  any  mesh  point  In  the  solution;  however,  when  approximating 
embeddec  shocks  as  Isentropic  discontinuities,  as 'done  in  CM3DT,  |DH|  may 
be  higher  than  0.05,  and  the  solution  may  still  be  judged  acceptable. 

For  solutions  that  were  judged  to  be  unacceptable,  it  was  noted 
that,  as  the  degree  of  the  indentation  increased,  the  flow  computed  along 
the  body  surface  frequently  reversed.  Any  calculation  with  this  flow 
reversal  at  the  surface  was  deemed  unacceptable.  (Since  the  calculations 
are  Inviscid,  this  flow  reversal  should  not  be  taken  as  an  indication  of 
separation,  which  Is  a  viscous  phenomenon.  It  is  probable  that  this  in- 
viscid  flow  reversal  is  indicative  of  the  limitations  of  the  approximation 
to  the  embedded  shock  used  in  the  CM3DT  analysis.) 
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Based  on  the  results  of  the  axisymmetric  calculations,  It  was 
determined  that  neither  freestream  Mach  number  nor  nose  radius  had  a 
significant  effect  on  the  ability  of  the  technique  to  compute  the  flow  on 
indented  shapes.  The  results  obtained  from  these  caluclaticns  could  most 
readily  be  correlated  In  terms  of  the  maximum  value  of  62  for  a  given  value 
of  0j  for  which  successful  calculations  could  be  made  (with  Rp  *  0).  A 
plot  of  92^x  versus  91  given  In  Figure  8.2. 

It  is  of  interest  to  note  that  the  curve  shown  In  Figure  8.2  can 
be  closely  approximated  by  a  straight  line  defined  as  02MAX  a  9l  +  15°. 

In  a  previous  study  evaluating  the  applicability  of  a  blunt  body  code 
(Reference  7),  it  was  determined  that  the  mesh  resolution  required  to  ob¬ 
tain  a  valid  finite  difference  solution  must  be  fine  enough  to  ensure  that 
no  mesh  interval  covers  a  region  of  flow  turning  angle  of  more  than  15°. 
Although  the  code  evaluated  In  Reference  7  is  applicable  only  to  convex 
geometries,  the  correspondence  of  the  A6  »  15°  limitation  arising  in  both 
the  CM3DT  evaluation  and  the  previous  effort  is  striking. 

Investigations  were  also  carried  out  to  assess  the  effects  of 
non-zero  fillet  radius  on  the  ability  of  the  technique  to  treat  indented 
nosetip  shapes.  In  general,  it  was  discovered  that  the  presence  of  a  fillet 
could  relax  the  02MAX  (Sj)  criterion  for  a  successful  calculation,  provided 
that  the  fillet  radius  is  large  enough  (and  the  mesh  resolution  fine  enough) 
to  have  mesh  points  in  the  finite  difference  grid  lying  on  the  fillet, 
located  such  that  A0  across  a  mesh  interval  is  no  larger  than  approximately 
15°. 


An  example  of  the  effects  that  can  be  expected  from  the  presence 
of  a  fillet  is  provided  by  the  Mildly  Indented  Body  (MIB),  an  axisymmetric 
indented  shape  described  in  Reference  5.  For  this  configuration,  with 
0i  =  26°  and  63  *  *6°,  the  criterion  shown  in  Figure  8.2  yields  ©2max  *  *10» 
with  a  zero  fillet  radius.  The  MIB  geometry,  however,  has  a  fillet  radius 
of  Rp  a  4.375  (with  R^  a  1.75),  and  the  CM3DT  code  is  known  to  give  good 
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Figure  8.2.  Limiis  of  Applicability 
the  CM3DT  Code  to  Indent 
Shapes 


results  for  this  case  (as  shown  in  Reference  1).  Thus,  the  presence  of 
the  fillet  has,  in  this  case,  extended  the  range  of  02  for  which  success¬ 
ful  calculations  can  be  expected. 

The  presence  of  a  sufficiently  large  fillet  cannot,  of  course, 
totally  relax  limitations  on  the  strength  of  the  compressions  that  can  be 
treated  (i.e.,  the  maximum  value  of  82  tor  which  a  calculation  is  feasible 
for  a  given  value  of  0i ).  In  general,  the  presence  of  a  fillet  should  be 
expected  to  Increase  the  maximum  allowable  value  of  @2  by  no  more  than  20%. 

The  Investigation  of  angle  of  attack  effects  on  axl symmetric  in¬ 
dented  nosetip  shapes  revealed  no  large  effects  relative  to  the  success 
of  CM3DT  calculations  for  a  given  geometric  configuration.  (It  was  noted, 
however,  that  as  the  angle  of  attack  Increased,  more  time  steps  were  re¬ 
quired  to  attain  a  satisfactory  solution.)  Because  of  the  almost  Infinite 
number  of  possible  geometric  variations,  no  attempt  was  made  to  study  angle 
of  attack  effects  on  asymmetric  indented  shapes. 

Based  on  the  results  of  this  effort,  the  applicability  of  the 
CM3DT  nosetip  code  to  a  given  Indented  nosetip  geometry  can  best  be  defined 
in  terms  of  the  versus  91  criterion  depicted  in  figure  8.2.  Use  of 

this  criterion  requires  the  user  to  approximate  a  particular  geometry  In 
terms  of  the  parameters  shown  in  Figure  8.1;  for  an  asymmetric  geometry, 
this  characterization  must  be  made  In  each  meridional  plane  to  be  computed. 
The  value  of  02^  c0rresP0n<i|in9  t0  the  particular  value  of  0j,  coupled 
with  the  potential  ameliorating  effects  of  a  non-zero  fillet  radius,  will 
provide  a  guide  to  the  user  as  to  whether  or  not  CM3DT  can  be  expected  to 
provide  a  successful  flow  field  calculation. 
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SECTION  10 


APPENDICES 


10.1 

GLOSSARY  OF  INPUT  VARIABLES 

Default  values  are  Indicated  In  parentheses  at  right. 

ACOF(I,J) 

1th  coefficient  In  jth  segment  of  analytic 
axlsynsnetrlc  body  geometry  definition 
(NSEG>0)  or  for  centerline  offset  In  pitch 
plane  (NSEG<0);  |NSEG|  segments 

ADIV(J) 

Abscissa  of  the  end  point  of  jth  segment 
and  the  beginning  point  of  tlie  jth  +  1st 
segment  of  analytic  axIsymmetrlcTody 
geometry  definition  (N$EG>0)  or  for 
centerline  offset  definition  (NSEG<0); 
j NSEG | -1  values  must  be  supplied 

AE(  I) 

Array  of  a  values  for  arbitrary  geometry 
option  with  centerline  offsets  (NPHI<0) 

ALPHA 

a,  Angle  of  attack  (degrees),  relative 
to  afterbody  centerline 

(0.0) 

ALT 

h.  Altitude  (feet).  If  IATMP  j*  0 

AMINF 

Mo,,  Freestream  Mach  number  (required 

If  QINF  not  Input) 

BE  ( I ) 

Array  of  b  values  for  arbitrary  geometry 
option  with  centerline  offsets 
(NPHI  »  -2,  -3,  -4) 

BETA 

B,  Sideslip  angle  (degrees),  relative 
to  afterbody  centerline 

(0.0) 
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CE(  I ) 

Array  of  c  values  for  arbitrary  geometry 
option  with  centerline  offsets 
(NPHI  ■  -3,  -4) 

DE(I) 

Array  of  d  values  for  arbitrary  geometry 
option  with  centerline  offsets  (NPHI  *  -4) 

DELA 

6ft,  Bend  angle  (degrees)  In  pitch  plane 
of  afterbody  axis 

(0.0) 

DELH 

Spacing  parameter  for  automatic  generation 
of  hinge  points 

(0.1) 

DELI 

6i«  Pitch  angle  of  orientation  of  nosetlp 
geometric  axis  relative  to  bent  after¬ 
body  axis  (degrees) 

(O.C) 

DEL2 

62,  Yaw  angle  of  orientation  of  nosetlp 
geometric  axis  relative  to  bent  after¬ 
body  axis  (degrees) 

(0.0) 

di(i) 

Array  of  centerline  offsets  In  pitch  plane, 
NPHI<0 

02(1) 

Array  of  centerline  offsets  in  yaw  plane, 
NPHI<0 

GIDEAL 

Y,  Ideal  gas  isentropic  exponent 

(1.4) 

GMW 

Gram  molecular  weight  of  equilibrium  air 

(28.9644) 

H1(I,L)* 

Location  of  1th  hinge  point,  x+iy  In  ith 
plane  for  IRTNGE  ■  1;  must  have  oC  * JA+2 
entries  In  each  plane,  NPHI  planes 

♦Complex  quantity 


IATMP 

Atmospheric  freestrcam  property  assignment 
Indicator: 

Freestream  pressure  and  temperature  from 

1962  standard  atmosphere 

(0) 

IHINGE 

Hinge  point  generation  indicator 
=0  Hinge  joints  generated  automatically 
=1  Hinge  points  input 

(0) 

I  LAM 

A-dlfferencing  calculation  indicator 
=0  for  non-conservation  calculation 
*1  for  A-scheme  calculation 

(0) 

IMAP 

Coordinate  transformation  indicator 
=0  For  ax i symmetric  mapping 
=1  For  asymmetric  mapping 

U) 

IRG 

Real  gas  option  indicator 

=0  For  ideal  gas  thermodynamics 

=1  For  equilibrium  air  thermodynamics 

(0) 

ISA 

Plane  of  symmetry  indicator 
=1  Pitch  plane  of  geometric  symmetry 
and  6=0 

-2  Asymmetric  shape  or  6^0 

(1) 

ITBOUT 

Shock  table  output  indicator  (IRG  =  1  only) 

=0  No  shock  table  output 
=1  Shock  tables  output 

(0) 

JA 

Number  of  "corners"  defined  by  hinge 
points  (maximum  7) 

JOUT 

Field  output  indicator 
=0  Field  data  only  output 
=1  Field  data  and  transformation  functions 
output 

=2  Field  data,  transformation  functions, 
and  circumferential  transformation  func¬ 
tions  output 

(0) 
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(10) 


KGRID 

KMAX 

KOUT 

KTAPE 

KSTART 

LMAX 

MMAX 

NMAX 

NPHI 


Number  of  time  steps  between  Inverse 
coordinate  transformations 


Maximum  number  of  time  steps 


Number  of  time  steps  between  complete  (50) 

output 


Number  of  time  steps  between  writing  output  (-1) 

on  tape 

<0  No  tape  output 


Step  number  for  restart  (-1) 

<0  Solution  initialized  at  step  0 


Number  of  circumferential  planes  in  (7) 

calculation  (maximum  9) 


Number  of  grid  points  along  body  In  each  (13) 

meridional  plane  (maximum  18) 


Number  of  grid  points  between  body  and  (6) 

shock  (maximum  11) 


Number  of  planes  of  arbitrary  body  geometry  (0) 

input  (maximum  9);  NPHI>0  for  no  offset 
in  geometric  centerline,  NPHI<0  for 
offset  geometric  centerline  option 
=1  Axlsymmetric  arbitrary  geometry 
-2  Circular  cross-sections  with  centers 
offset  from  geometric  centerline,  pitch 
plane  of  geometric  symmetry  assumed 
=3  Bi-elliptic  cross-sections,  pitch  plane 
of  geometric  symmetry  assumed 
=4  and  ISA  =  2  Quad-elliptic  cross-sections 
=4  and  ISA  *  1,  or  NPHI >4  Arbitrary  cross-sections 
=-l  Circular  cross-sections  with  offset 
geometric  centerline 
=-2  Elliptic  cross-sections  with  offset 
geometric  centerline 

*-3  Bl-elliptlc  cross-sections  with  offset 
geometric  centerline 

*-4  Quad-elliptic  cross-sections  with  offset 
geometric  centerline 
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NPTS(L) 


t  v 5 v-?*¥S- ifi "V.’ v/ /v-ff ; V’’*“75 ' 


NPHI>0,  number  of  (XB,YB)  points  Input  in 
PHIBD(L)  plane  for  arbitrary  geometry 
option  with  no  centerline  offset 
(maximum  value  of  50) 

NPHI<0,  number  of  entries  in  XE,  AE,  BE, 

CE,  DE,  Dl,  and  D2  arrays  for  arbitrary 
geometry  option  with  centerline  offset; 
only  NPTS(l)  required  (maximum  value  of  50) 


NSEG  >0  Number  of  segments  used  to  define  (0) 

analytic  axi symmetric  body  geometry 
(maximum  6) 

<0  For  NPHI<0,  number  of  segments  used 
to  define  pitch  plane  centerline  offsets 
for  arbitrary  body  geometry  (maximum 
value  of  | NSEG |  is  6) 


NUWTB  Numoer  of  entries  in  shock  table  (91) 

(real  gas) 


PINF  p„,  Freestream  static  pressure  (real  gas,  (1.0) 

lbf/ft2) 


PHIBD(L)  Circumferential  planes  for  arbitrary  geometry 

input  (degrees);  NPHI>0  (maximum  of  9) 


QINF  Voo,  Freestream  velocity  (required  If  AMINF  (0.0) 

not  input)  (real  gas,  ft/sec) 

RIDEAL  R,  Gas  constant  (ideal  gas)  (1.0) 

RN  Rn,  Effective  nose  radius,  used  to  estimate  (0.5) 

Initial  shock  stand-off  distance  and  locate 
hinge  points,  (IHINGE  =  0) 

STAB  Stability  factor  in  step-size  calculation  (1.0) 

TEST(l)  Convergence  criterion  for  geometry  iteration  (10~6) 

TEST(3)  e  for  hinge  points  on  horizontal  axis  (10-6) 
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TEST ( 8 ) 


Factor  on  shock  stand-off  distance  for  (1,5) 

automatic  generation  of  H1(1,L) 


TEST(ll)  Convergence  criterion  for  shock  property  (10"') 

iteration  (real  gas) 

TEST (12)  Convergence  criterion  for  stagnation  (10~6) 

conditions  iteration  (real  gas) 

TINF  Too,  Freestream  static  temperature  (1.0) 

(real  gas,  °R) 

UWTBMN  Minimum  value  of  u/V^**,  independent  (0.1) 

variable  for  shock  tables 

UWTBMX  Maximum  value  of  u/V«**,  Independent  (1.2) 

variable  for  shock  tables 


XB(I,L)  Abscissa  of  1th  point  in  arbitrary  geometry 

description  in  plane  PHIBD(L ) ;  NPHI>0 
(XB( 1  ,L)  *  XBO) 


XBB(L)  Abscissa  of  body  point  in  plane  PHIBD(L) 

defining  downstream  boundary  of  the 
computational  region  (NPKI  values  to  be 
specified  if  NPHI>0,  one  value  if  NPHI<0 
or  if  NSEG>0) 


XBO  Location  of  body  point  on  nosetip  (0.0) 

geometric  axis 


XE ( I )  Array  of  9  values  for  arbitrary  geometry 

option  with  centerline  offsets  (NPHI<0) 
(NPTS(l)  values) 


XHO  Value  of  9  where  nosetip  geometric  axis 

intersects  bent  afterbody  axis 


**u  is  the  normal  component  of  the  freestream  velocity  relative  to 
the  shock 


YB(I.L) 


Ordinate  of  ith  point  In  arbitrary  geometry 
description-!)!  plane  PHIBD(L);  NPHI>0 
(YB(1,L)  *  0.0) 


Z I DEAL 


z.  Compressibility  factor  (Ideal  gas) 


GLOSSARY  OF  OUTPUT  VARIABLES 


10.2 

AG 

AL(L) 

AMSH 

CIS1 

CIS2 

DELTA 

DH 

DYDPH 

DYDX 

ETA 

ETAB 

ETABTH 

ETABXI 

ETAS 

ETAST 


1 9 1  *  1 35/3*1 

a(<j>).  Stretching  parameter 

Mach  number  downstream  of  normal  shock 

C  »  cos(-o)) 

S  *  sin(-<u) 

Shock  layer  thickness  at  the  centerline 

(H-Ho)/Hq,  Variation  of  total  enthalpy 

3y./34>,  Circumferential  body  surface  derivative 
In  physical  space 

3yb/3x,  Longitudinal  body  surface  derivative 
In  physical  space 

n»  Ordinate  in  transformed  space 

b(E»9)»  Body  surface  in  transformed  space 

3b/38,  Circumferential  body  surface  derivative 
in  transformed  space 

3b/ 3E.  Longitudinal  body  surface  derivative 
In  transformed  space 

c(E,6,t),  Bow  shock  surface  In  transformed  space 
3c/3t,  Shock  velocity  in  transformed  space 
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ETASTH 

ETASXI 

ETSOGM 

ETSTOG 

G 

H 

H(I,J,L) 

HTOT 

IMG 

IMGPHI 

IMZETPHI 

M 

MACH 

P 

PHI 


3c/38,  Circumferential  shock  surface  derivative 
In  transformed  space 

3c/3C.  Longitudinal  shock  surface  derivative 
in  transformed  space 

Maximum  absolute  value  of  ETSTOG  at  all  shock 
points  at  a  giver,  time  step 

|^-|(3c/3t),  Shock  velocity  In  physical  space 

y,  Isentroplc  exponent 

h.  Static  enthalpy 

hi,j,k>  Image  of  1th  hinge  point  after  (j-l)th^ 
transformation  Tn  Jtth^  41  plane 

Hq,  Frees tream  total  enthalpy 

Im(g) 

ImV 

In,(^) 

M  ,  Freestream  Mach  number 
00' 

M,  Local  Mach  number 
p,  Static  pressure 

4>,  Circumferential  coordinate  in  physical  space 
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.  „„U  vtt'v  „„  .  _  _  . . 


Re($) 


Im(<t) 


POWER(J.L)  5j,i»  Exponent  of  jth  transformation  In  Ath 
$  plane 


Static  pressure  downstream  of  normal  shock 
(Real  gas,  lbf/ft2) 


PSTAG 


p_,  Stagnation  pressure  (Real  gas,  lbf/ft2) 


V^,  Freestream  velocity  (Real  gas,  ft/sec) 

^PooVoo2,  Freestream  dynamic  pressure  (Real  gas, 
*  lbf/ft2) 


Re(g) 


REGPHI 


Re<V 


REZETPHI 


ReV 


p,  Density  (Real  gas,  slugs/ft3) 


s/R,  Non-dimensional  entropy 


Non-dimensional  entropy  behind  normal  shock 


T,  Static  temperature  (Real  gas,  °R) 


THETA 


8,  Circumferential  coordinate  in  transformed  space 


t.  Time  coordinate 


u 

V 

w 

WRMS 

X 

XI 
XIB 

XIL(L) 

XILTH(L) 

XIS 

XB 

Y 

YB 

YS 


u,  5-veloclty  component 

v,  n-veloclty  component 

w,  9 -velocity  component 

Root -mean-square  of  shock  velocities  In  physical 
space 

x,  Axial  coordinate  In  physical  space 
5,  Abscissa  in  transformed  space 

5  at  body  point 

5L(e),  Downstream  boundary  of  computational  space 
d5L(Q)/d6 

5  at  bow  shock  point 
x  at  body  point 

y,  Radial  coordinate  in  physical  space 
y  at  body  point 

y  at  bow  shock  point 


'tfW*  JWjt- 
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GLOSSARY  OF  SUBROUTINE  NAMES 


Routines  indicated  by  (*)  were  obtained  from  the  flow  field 
code  described  in  References  2  and  3. 


AIR* 

ATERR 

ATMP* 


BFL3CM 

BFR3CM 

BFW3CM 

CM3DT 

CM3ST 


CUFT1* 


CUFT2* 

FREE* 


GE0M3 

GE0M4 


GE0M4B 


Tables  of  equilibrium  air  thermodynamic  properties 


Prints  out  error  message 


Determi nati on  of  freestream  pressure  and 
temperature  from  tables  of  1962  Standard  Atmosphere 


Binary  file  routine 


Reads  binary  file  (entry  point  in  BFL3CM) 


Writes  binary  file  (entry  pbint  in  BFL3CM) 


Main  routine;  initialization 


Calculates  starting  line  data  fo1"  afterbody 
computation 


Cubic  spline  fit  routines 


Calculates  freestream  conditions 


Analytic  body  geometry  subroutine 


Initialization  for  arbitrary  body  geometry 
definition,  with  no  centerline  offset 


Arbitrary  body  geometry  routine  with  no 
centerline  offset  (entry  point  in  GE0M4) 
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GE0M5 

GE0M5B 

GRID3D 

HING3D 

IDEAL* 

LINB* 

MAP3D 

0UT30 

PVE3DNC 

PVE3DL 

REALP 

RGAS+ 


Initialisation  for  arbitrary  body  geometry 
definition,  with  centerline  offset 


Arbitrary  body  geometry  routine  with 
centerline  offset  (entry  point  in 
GE0M5) 


Calculates  coordinate  transformation  para¬ 
meters  at  grid  points 


Calculates  hinge  point  locations  and 
transformations  of  hinge  points 


Calculates  ideal  gas  properties 

Solves  system  of  linear  equations 
(used  in  spHne  fit) 


Conformal  mapping  transformation  function 
(complex) 


Output  routine 


Time -dependent  calculations  in  non-conservation 
form 


Time-dependent  calculations  with  the  ^.-differencing 
scheme 


Determines  real  part  of  a  complex  number 
(used  in  CM3 ST) 


Calculates  real  gas  properties 


RIERR  Prints  out  error  message  if  range  of  real  gas 

tables  is  exceeded 


+  This  routine  Is  the  same  as  the  REAL  routine  in  Reference  3:  the  name  has 
been  changed  to  allow  use  of  the  Intrinsic  function  REAL  in  complex  arith¬ 
metic  with  CM3DT. 
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SHKTAB*  Constructs  shock  tables 

SHOCK*  Calculates  properties  downstream  of  shock  for 

specified  nominal  frees tream  velocity  component 

STAGPT*  Calculates  stagnation  point  conditions 

TBL*  Table  search 

TBL1*  Single  table  look-up 

TLU1*  Table  look-up  (used  In  spline  fit) 
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10.4 


SUMMARY  OF  COORDINATE  TRANSFORMATIONS 


To  provide  maximum  flexibility  to  the  user,  the  CM3DT  nosetlp 
flow  field  code  has  been  formulated  to  allow  definition  of  the  nosetlp 
geometry  In  a  coordinate  system  that  Is  convenient  to  the  user  and  Is  not 
necessarily  aligned  with  the  vehicle  afterbody  coordinate  system.  In 
addition,  provision  has  been  made  for  a  bent  afterbody  axis  to  further 
generalize  the  applicability  of  the  CM3DT  code.  This  section  details 
the  relationships  between  the  various  coordinate  systems  that  are  used  In 
the  CM3DT  code. 

The  afterbody  cylindrical  coordinate  system,  (r,6,z),  Is  defined 
as  shown  in  Figure  10.1,  where  the  axis  may  be  bent  In  the  pitch  plane  by 
an  angle  6^  at  z  *  z.  Defining  an  afterbody  Cartesian  frame  (x^y^z^) 
oriented  with  the  unbent  afterbody  axis,  the  Cartesian  coordinates  may  be 
expressed  as 

Xj  *  z  -  7  (10.1) 

yj  *  r  cos  9  +  d(z)  (10.2) 

Zj  «  r  sin  9  (10.3) 

where 


d(z)  *  (z  -  z)  tan  6^  .  (10.4) 

Another  Cartesian  coordinate  system  (x'.y'.z1)  may  be  defined 
that  is  oriented  with  the  bent  afterbody  axis.  This  system  is  related  to 
the  (xj.y^.Zj)  frame  by 


X  1  ' 

r-  COS  6^ 

-  sin  6a 

0 

'  V 

s  y*  -  * 

sin  «A 

cos  6 

0  : 

'  yl  ; 

L  z‘  i 

0 

0 

1  - 

1  *1- 
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The  next  coordinate  system  of  interest  is  the  nosetip  geometric 
cylindrical  system  (x,y,$),  in  which  the  nosetip  geometry  Is  to  be  defined 
in  CM3DT.  The  associated  Cartesian  system,  denoted  by  (xc,yc,zc),  is 
defined  from 


A 

-  X 

(10.6) 

=  y  cos  $ 

(10.7) 

=  y  sin  3 

(10.8) 

As  discussed  in  Section  4  of  this  report,  the  nosetip  geometric 
axis  intersects  the  bent  afterbody  axis  at  z  =  zQ  and  x  =  xQ,  and  may  be 
rotated  in  both  pitch  and  yaw  by  the  angles  fij  and  62,  respectively,  as 
shown  in  Figure  10.2.  The  (xr,yr,z_)  and  (x'.y'.z1)  Cartesian  systems  may 
then  be  related  through 


A  A 

Vx0l  ! 

1  i 

rc.os  61  -  sin  Sj 

0-, 

“COj  8^ 

0 

-sin 

'X'  +  (z-zQ)/cos  fiA'| 

?e  }■ 

sin  <5.  cos  <5. 

I  1  1 

0 

0 

1 

0 

y 

Sc  J‘ 

o 

o 

— 1 
rH 

-sin  &2 

0 

cos  6^ 

(10.9) 


The  (x,y,$)  nosetip  geometric  coordinate  system  also  serves  as 
the  ( x ,y ,cp )  computational  coordinate  system  in  the  CM3DT  code  for  those 
cases  where  there  is  no  offset  to  the  cross-section  centers  (NSEG>0  or 
NPHI>0).  However,  if  the  cross-section  centers  are  offset  (NPHI<0 )  the 
computational  axis  will  not  be  coincident  with  the  geometric  axis,  and 
additional  transformations  will  be  required. 
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The  nosetip  computational  axis  is  defined  as  a  straight  line 
that  passes  through  the  points  (xbo,  d1(J,  d20),  (xQ,  0,  0)  in  the  geometric 
Cartesian  coordinate  system,  where  xbo  *  XBO,  diQ  =  Dl(l),  d20  =  02(1), 
and  xQ  =  XHO.  As  shown  in  Figure  10.3,  the  orientation  of  the  computa¬ 
tional  axis  to  the  geometric  axis  can  be  described  by  the  angles  and 

62n: 

<$2N  =  d2o/(x0  -  xb())  (10.10 

tan  61n  =  dlo  cos  52N/(*0  -  xbQ)  .  (10.11 

Defining  the  Cartesian  frame  oriented  with  the  computational  axis 
as  (xc,yc,zc),  the  transformation  to  the  computational  frame  becomes 

xc~(V^bo)/(cos  61n  cos  62n) *1  r  cos  61n  -sin  61n  0- 

yc  *  sin  fiiN  cos  62n  o  : 

z.  j  L  0  0  1  - 

-  cos  62n  0  -sin  52n  -j  f  xc-xQ  y 

0  1  0  j  y.  '  .  (10.12 

j  Ac  [ 

_  sin  62n  0  cos  62n  1  zc  j 

This  transformation  is  defined  such  that  the  most  forward  body  point  in 

thq  geometric  system,  x  ®  xb(),  has  the  value  xc  »  0  in  the  computational 
system.  The  CM3DT  code  automatically  redefines  the  input  parameters 
XBO  and  XHO  to  be  0  and  (x  -  xbo)/cos  6^  cos  62^,  respectively,  when 
a  new  computational  axis  Is  defined.  In  addition,  the  code  also  redefines 
the  angles  <5j  and  62  so  that  they  represent  the  orientation  of  the  computa¬ 
tional  axis  to  the  bent  afterbody  axis. 


The  redefinition  of  6^  and  62  is  accomplished  by  expressing 
the  coordinates  of  the  point  (xbQ,  dj0,  d20)  in  the  geometric  coordinate 
system  in  terms  of  the  Cartesian  system  oriented  with  the  bent  afterbody 
axis  (by  inversion  of  Equation  (10.9))  to  yield 

0  sin  -> 

1  0  x 

0  cos 

..  COS  sin  61  0  ,  Xb0  -  xQ  , 

.  -sin  cos  0  dlo  |  .  (10.13) 

i  » 

L  0  0  1  J  L  d2o 

(Note  that  in  Equation  (10.13)  the  angles  and  52  are  the  values  input 
by  the  user  to  define  the  orientation  of  the  geometric  axis.) 

New  values  of  6j  and  <$2  (defining  the  orientation  of  the  compu¬ 
tational  axis)  can  be  determined  from  the  results  of  Equation  (10.13)  as 

tan  <$2  =  -z ' /{x ’  -  ( z  -zQ)/cos  5ft}  (10.14) 

tan  61  =  -y'  cos  <52/{x'  -  (z  -  z0)/cos  6 A)  .  (10.15) 
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SAMPLE  CM30T  OUTPUT 


SAMPLE  CASE  1.  TITLE  SHEET. 

SAMPLE  CASE  2.  TITLE  SHEET. 

SAMPLE  CASE  3.  TITLE  SHEET. 

SAMPLE  CASE  3  (AFTERBODY).  TITLE  SHEET. 

ANALYTIC  GEOMETRY  INPUT  (SAMPLE  CASE  1). 

ARBITRARY  GEOMETRY  INPUT  (NPHI  *  1)  (SAMPLE  CASE  2). 
ARBITRARY  GEOMETRY  INPUT  (NPHI  =  2)  (SAMPLE  CASE  3). 
HINGE  POINT  TRANSFORMATION  (SAMPLE  CASE  3). 
CIRCUMFERENTIAL  TRANSFORMATION  DATA  (SAMPLE  CASE  3). 
BODY  GEOMETRY  INFORMATION  (SAMPLE  CASE  2). 

SAMPLE  TIME  HISTORY  OUTPUT  (SAMPLE  CASE  1). 

FIELD  OUTPUT  (SAMPLE  CASE  l). 

FIELD  OUTPUT  (SAMPLE  CASE  2), 

FIELD  OUTPUT  (SAMPLE  CASE  3). 

FIELD  OUTPUT,  JOUT  =  1  (SAMPLE  CASE  3). 

FIELD  OUTPUT,  JOUT  =  2  (SAMPLE  CASE  3). 
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F«F £  STREAK  CONDITIONS 


HINCE  POINTS  CFNEMTtO  AO  10HAI ICALIT 
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SAMPLE  CASE  3.  TITLE  SHEET 


EMFRHAL  PARAHE TERS-KENSF- 


d.)  SAMPLE  CASE  3  (AFTERBODY).  TITLE  SHEET 
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ANALYTIC  GEOMETRY  INPUT  (SAMPLE  CASE  1). 


f.)  ARBITRARY  GEOMETRY  INPUT  (NPHI  =  1)  (SAMPLE  CASE  2). 
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.(cont'd).  FIELD  OUTPUT,  JOUT  =  1  (SAMPLE  CASE  3  ). 
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